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It  tUPPUNINTARV  NOTie 


A  general  circulation  model  for  climate,  emphasizing  physical  phenomenology,  ii  neve  loped  as  p<*rt  of  an 
overall  consideration:  whether  it  is  feasible  to  control  or  advertently  modify  the  earth’s  surface  climate  at 
the  middle  and  high  northern  latitudes.  Ultimately,  the  answer  to  this  feasibility  question  will  be  used  tc 
evaluate  theoretically  the  current  capabilities  of  world  superpowers  in  their  possible"  utilisation  of  controlled 
climate  modification  for  political  or  military  purposes. 

The  use  of  the  model,  along  with  the  most  recent  reduction  of  experimental  data  for  the  long-term  statistics 
of  the  current  general  circulation,  allows  the  semi -experimental  determination  of  certain  statistics  which  are 
difficult  to  measure  directly.  These  are  the  wave  and  amplitude  parameters  of  the  baroclinic  disturbances, 
and  the  vertical  transports,  which  the  model  determines  as  a  function  of  latitude.  A  local  test  of  the  model 
proved  successful:  for  the  latitude  where  the  large-scale  transport  of  angular  momentum  vanishes  as  the  large- 
scale  transport  of  sensible  heat  reaches  a  maximum,  experimental  and  semi-experimental  results  show  agreement. 

The  indicates  that  the  transports  of  angular  moment-im,  sensible  heat,  and  specific  humidity  do  depend 

on  the  zonal  (eastward  direction)  circulation.  This  resolves  the  puzzle  of  previous  results  in  which  wave  and 
amplitude  parameters  of  the  baroclinic  disturbances  were  assumed  to  be  constant.  Also,  the  model  indicates 
that  these  transports  are  coupled  to  the  vertical  ones.  Only  one  eigenvalue  is  required  by  the  model.  This 
involves  the  static  stability  of  the  troposphere. 
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FOREWORD 


This  paper  considers  the  question  of  whether  controlled  or  ad¬ 
vertent  modification  of  the  earth's  climate  at  the  middle  and  high 
northern  latitudes  is  feasible.  The  work  was  performed  by  IDA  for 
the  Defense  Advanced  Research  Projects  Agency  (DARPA). 

The  author  is  indebted  to  Dr.  W.L.  Gates  of  the  Rand  Corporation 
for  his  comprehensive  review  of  the  manuscript  of  this  paper,  and  to 
Dr.  W.M.  Washington  of  the  National  Center  for  Atmospheric  Research 
for  important  general  comments  on  the  main  subject  of  this  paper. 

One  of  Dr.  Washington's  worthwhile  comments  concerns  the  development 
of  numerical  models  of  the  atmosphere  for  the  prediction  of  climate 
changes,  i.e.,  "it  is  the  feeling  among  most  modellers  that  neither 
simplified  nor  complex  models  are  developed  to  the  point  of  giving 
definite  answers  to  climate  change." 

The  author  acknowledges  the  support  of  his  colleagues  at  IDA, 
i.e.,  Dr.  Ernest  Bauer  for  reviewing  the  manuscript  and  Mr.  Dennis  F. 
Deriggi  for  his  work  on  the  numerical  calculations  for  this  paper. 

The  author  also  thanks  Dr.  Robert  Vaglio-Laurin  for  some  valuable 
technical  discussions  concerning  two-dimensional  turbulence. 
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SUMMARY 


This  paper  considers v the  feasibility  of  controlled  or  advertent 
modification  of  the  climate  of  the  earth’s  surface  at  the  middle  and 
high  northern  latitudes.*  The  ultimate  aim  of  this  feasibility  ques¬ 
tion  is  the  theoretical  evaluation  of  the  current  capabilities  of  a 
world  superpower  to  use  controlled  climate  modification  for  political 
or  military  purposes.  The  timeliness  of  this  feasibility  question 
stems  from  two  technological  developments  which  will  enhance  the  un¬ 
derstanding  of  tropospheric  phenomena:  (1)  the  advent  of  meteorolog¬ 
ical  satellites,  which  provides  a  capability  for  worldwide,  near  real¬ 
time  measurements  of  the  main  characteristic  parameters  of  the  tropo¬ 
sphere  (e.g.,  the  temperature  field,  humidity  etc.);  and  (2)  the  de¬ 
velopment  of  large-memory,  high-speed  electronic  computers,  which  will 
allow  the  reduction  of  large  amounts  of  data  and  the  performance  of 
theoretical  calculations  with  increasing  space  resolutions  both  over 
the  earth’s  surface  and  along  the  vertical  dimension  of  the  tropo¬ 
sphere. 

Considerations  of  the  feasibility  of  advertent  climate  modifica¬ 
tion  at  the  middle  and  high  northern  latitudes  must  strive  for  a  maxi¬ 
mum  understanding  of  the  basic  processes  that  control  atmospheric  and 
oceanic  phenomena  near  the  earth's  surface.  The  study  of  such  proc¬ 
esses  during  recent  years  has  emphasized  the  development  of  complex 
numerical  models  of  the  troposphere  using  increasing  space  resolutions; 
a  trend  that  requires  the  use  of  large  electronic  computers  of  ever 


The  concept  of  controlled  or  advertent  climate  modification  is  dif¬ 
ferent  from  that  of  either  inadvertent  climate  modification  caused 
by  the  impact  of  modern  man*s  activities  on  his  environment  or  ad¬ 
vertent  weather  modification. 


v 
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increasing  capabilities.  Such  procedures  yield  a  wealth  of  numerical 
data  which  do  not  necessarily  enhance  the  understanding  of  the  con¬ 
trolling  processes.  The  overall  considerations  of  the  feasib'' lity  of 
advertent  climate  modification  have  therefore  used  two  approaches: 

(1)  an  attempt  to  capitalize  on  the  availability  of  the  complex  numer¬ 
ical  models  for  the  general  circulation  of  air  within  the  troposphere 
(e.g.,  the  Mintz-Arakawa  model,  etc.),  and  (2)  an  investment  of  a 
limited  manpower  level  in  developing  simpler  complementary  models  of 
the  general  circulation  for  climate  that  emphasize  physical  phenome¬ 
nology.  This  report  is  part  of  this  second  approach. 

The  fundamental  question  concerning  the  feasibility  of  advertent 
climate  modification  is  broken  down  in  this  report  into  four  basic 
subquestions : 

1.  What  are  the  key  processes  that  control  the  climate  at  the 
earth’s  surface? 

2.  Can  the  results  of  these  processes  be  described  uniquely  from 
the  use  of  the  physical  laws  controlling  such  processes? 

3.  What  is  the  sensitivity  of  the  climate  at  the  earth’s  surface 
to  changes  in  these  key  processes? 

4.  Does  man  have  the  technology  now  to  use  such  a  sensitivity  to 
control  his  climate? 

A  positive  answer  to  the  fundamental  question  of  the  feasibility 
of  advertent  climate  modification  requires  conclusive  positive  answers 
t°  each  of  these  subquestions.  Thus,  the  answer  to  the  first  subques¬ 
tion  must  take  the  form  of  an  accurate  model  for  climate,  i.e.,  a 
mathematical  description  of  the  basic  processes  that  control  the  long¬ 
term  (e.g.,  annual,  seasonal,  monthly)  statistics  of  the  general  cir¬ 
culation.  The  second  subquestion  is  the  most  difficult  to  answer  con¬ 
clusively,  for  it  requires  as  yet  unknown  ergodic  criteria  for  the 
ruling  out  of  all  but  one  of  any  multiple  solutions  from  the  model 
that  predicts  changes  in  the  earth’s  climate.  The  answer  to  the  third 
subquestion  must  be  based  on  the  previous  ones,  i.e.,  on  the  use  of  an 
accurate  model  that  yields  unique  predictions  for  changes  in  the  earth’s 
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climate;  the  earth’s  climate  must  then  have  an  appropriate  degree  of 
sensitivity  to  changes  in  the  key  processes.  The  answer  to  the  final 
subquestion  must  be  based  on  the  answer  to  the  third  subquestion. 
Because  of  the  time  lag  between  initiation  and  completion  of  a  major 
project,  which  is  usually  a  major  fraction  of  a  decade,  a  positive 
answer  to  the  fourth  subquestion  becomes  necessary  if  man  is  to  con¬ 
trol  his  climate  as  early  as  in  the  1980s. 

The  results  of  this  paper  are  addressed  mainly  to  the  first  sub¬ 
question,  i.e.,  the  development  of  a  model  of  the  general  circulation 
appropriate  for  climate  by  emphasizing  physical  phenomenology.  For 
this  reason,  the  scope  of  this  report  attempts  (1)  to  present  a  rather 
unified  picture  of  the  main  developments  in  meteorology  during  the 
last  30  years,  which  have  enhanced  the  understanding  of  the  basic 
processes  that  control  the  earth's  climate,  and  (2)  to  describe  the 
physical  and  mathematical  arguments  in  the  development  of  a  model, 
based  on  physical  phenomenology,  with  enough  detail  to  bring  out  as 
clearly  as  possible  all  the  assumptions  used  in  such  a  model. 

Key  processes  that  control  the  earth’s  climate,  as  determined  by 
the  long-term  statistics  of  the  general  circulation  of  air  within  the 
troposphere,  are  the  poleward  transports  of  angular  momentum  and  sen¬ 
sible  heat  (or  enthalpy).*  These  processes  involve  two  different  but 
coupled  types  of  phenomena: 

1.  Net  heating  rates  of  the  troposphere  as  a  function  of  lati¬ 
tude  by  solar  and  terrestrial  radiation,  condensation  of 
water  vapor,  and  eddy  fluxes  of  sensible  and  latent  heat  to 
the  troposphere  from  the  surface  boundary  layer  (about  1  km 
thick)  next  to  the  earth's  surface  (e.g.,  Figs.  1  and  2, 
Section  III). 

2.  Large-scale  turbulent  or  eddy  motion  that  is  triggered  by  the 
net  heating  rates  of  the  troposphere;  because  of  the  condi¬ 
tion  of  very  small  vertical  winds  or  hydrostatic  equilibrium 

* 

The  troposphere  is  the  lowest  layer  of  the  atmosphere,  about  16  km 
thick  at  the  equator  and  about  8  km  at  the  poles. 
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in  the  direction  vertical  to  the  earth’s  surface,  the  tur¬ 
bulent  motion  is  two-dimensional  with  its  own  characteristic 
spectra  (Fig.  4,  Section  ill).  However,  in  contrast  with 
the  problem  of  long-range  weather  forecasting,  only  low 
hemispherical  wave  numbers  (i.e.,  number  of  wavelengths  in  a 
latitudinal  circle)  are  of  interest  for  climate. 

The  eddy  motion  at  low  wave  numbers  produces  the  large-scale 
poleward  transports  of  angular  momentum  and  sensible  heat,  which  in 
turn  set  up  a  meridional  cell  structure  of  a  statistical  nature  for 
the  general  circulation  (e.g.,  Fig.  5,  Section  III).  The  number  of 
the  meridional  cells  in  each  hemisphere  varies  between  two  and  three, 
which  corresponds,  respectively,  to  the  two  extremes  of  interglacial 
and  glacial  climates  in  the  current  Pleistocene  Epoch.  Large-scale 
transports  of  angular  momentum  and  sensible  heat  take  place  through 
a  mechanism  that  involves  the  presence  of  streamlines  with  an  asym¬ 
metrical  sinusoidal  form,  which  are  characterized  by  the  associated 
NE-SW  tilt  (a)  from  the  meridian  of  their  trough  line  (Fig.  12,  Sec¬ 
tion  VII).  This  mechanism  produces  transports  of  angular  momentum 
that  are  countergradient  at  middle  and  higher  latitudes,  i.e.,  trans- ' 
port  that  cannot  be  expressed  by  the  use  of  an  eddy  viscosity  and  the 
latitudinal  gradient  of  the  average  flow.  Since  the  transport  of 
angular  momentum  is  largest  in  the  upper  half  of  the  troposphere,  a 
model  of  the  general  circulation  for  climate  can  be  set  up  by  using 
a  two- level  model  of  the  troposphere  (Fig.  13,  Section  VII);  i.e., 
the  flow  pattern  in  the  upper  half  of  the  troposphere  may  be  repre¬ 
sented  by  that  at  the  250-mb  level,  while  the  flow  pattern  in  the 
lower  half  of  the  troposphere  is  represented  by  that  at  the  750-mb 
level.  The  conservation  equations  are  applied  at  these  levels  together 
with  the  analytical  formulation  of  the  flow  ( if* )  and  temperature  (^*) 
disturbances  induced  by  the  eddy  motion  (Fig.  12).  Based  on  the  Me¬ 
ridional  cell  structure  of  the  general  circulation  for  the  two  extreme 
conditions  of  the  interglacial  and  glacial  climates  of  the  current 
Pleistocene  Epoch,  a  further  assumption  used  in  the  model  is  that  the 
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role  of  the  angular  momentum  transport  is  always  more  important  in  the 
upper  half  as  compared  with  that  in  the  lower  half  of  the  troposphere. 

Specific  results  concerning  the  first  subquestion  are  then  as 
follows: 

1.  A  ruodel  of  the  general  circulation  for  advertent  climate 
modification  has  been  formulated  based  on  a  variable  tilt 
from  the  meridian  for  both  the  flow  and  temperature  distur¬ 
bances,  a  variable  phase  lag  between  the  flow  and  tempera¬ 
ture  disturbances,  and  a  variable  amplitude  of  these  non¬ 
linear  baroclinic  disturbances. 

2.  The  use  of  this  model  based  on  variable  wave  and  amplitude 
parameters  of  the  baroclinic  disturbances,  together  with  the 
most  recent  reduction  of  experimental  data  for  the  long-term 
statistics  of  the  current  general  circulation,  allows  the 
semi -experimental  determination  as  a  function  of  latitude  of 
statistics  that  are  difficult  to  measure  directly,  i.e.,  the 
wave  and  amplitude  parameters  of  the  baroclinic  disturbances 
as  well  as  the  vertical  transports. 

3.  A  successful  local  t^est  of  the  fundamental  formulation  of  the 
nonlinear  baroclinic  disturbances  is  obtained  from  the  agree¬ 
ment  between  experimental  and  semi -experimental  results  for 
the  latitude  where  the  large-scale  transport  of  angular  mo¬ 
mentum,  vanishes  as  the  large-scale  transport  of  sensible  heat 
reaches  a  maximum  value. 

4.  The  model  based  on  variable  wave  and  amplitude  parameters  of 
the  baroclinic  disturbances  indicates  that  the  transports  of 
angular  momentum,  sensible  heat,  and  specific  humidity  do  de¬ 
pend  on  the  zonal  (i.e.,  eastward  direction)  circulation,  a 
result  that  resolves  the  puzzle  of  previous  results  assuming 
constant  wave  ard  amplitude  parameters  of  the  baroclinic  dis¬ 
turbances.  The  model  based  on  variable  wave  and  amplitude 
parameters  further  indicates  that  the  foregoing  transports 
are  also  coupled  to  the  vertical  transports. 
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5.  The  model  based  on  variable  wave  and  amplitude  parameters  of 
the  baroclinic  disturbances  requires  only  one  eigenvalue 
parameter,  which  involves  the  static  stability  of  the  tropo¬ 
sphere.  This  result  implies  that  there  is  no  need  to  use-- 
as  in  the  model  using  constant  amplitude  and  wave  parameters- 
a  weak  boundary  condition  at  the  north  pole  to  establish  the 
magnitude  of  the  transports  of  angular  momentum  and  sensible 
heat  at  latitudes  lower  than  those  near  the  north  pole. 

Integrations  of  this  model  can  yield  the  changes  in  the  basic 
statistics  that  control  the  general  circulation  within  the  troposphere 
as  a  function  of  essentially  the  advertent  modifications  of  the  pole- 
ward  heat  flux  at  the  middle  and  high  latitudes;  a  modification  that 
would  result  from  any  advertent  change  of  the  surface  albedo,  for  ex¬ 
ample,  and  its  related  impact  on  the  net  heating  rates  of  the  tropo¬ 
sphere  at  these  latitudes.  These  integrations  can  use  known  boundary 
conditions  at  the  northern  edge  of  the  unperturbed  meridional  cell  at 
low  latitudes,  i.e.,  the  Hadley  ceU,  in  which  the  ITCZ  phenomena  play 
a  dominant  role  (Fig.  7,  Section  III).  The  changes  in  the  long-term 
statistics  determinable  from  the  model  involve  the  large-scale  pole- 
ward  transport  of  angular  momentum  at  the  250-mb  level,  the  mean  sen¬ 
sible  heat  and  specific  humidity  transports,  the  vertical  transports, 
and  the  meridional  as  well  as  zonal  eddy  kinetic  energies  in  the  mid¬ 
dle  and  in  the  upper  half  of  the  troposphere. 

The  second  basic  subquestion  identified  earlier  has  to  do  with 
the  uniqueness  (or  lack  of  it)  of  a  predicted  change  of  the  climate  at 
the  earth’s  surface  as  derived  from  the  nonlinear  conservation  equa¬ 
tions.  This  problem  is  again  different  from  the  predictability  of 
weather  forecasting  for  a  few  days,  since  the  fundamental  question 
here  concerns  the  possible  multiplicity  of  solutions  for  given  non¬ 
linear  differential  equations  with  appropriate  boundary  conditions. 

If  such  solution  is  unique,  then  there  is  only  one  prediction  for  the 
change  of  the  climate;  the  solution  then  becomes  deterministic  or 
transitive.  If  there  are  multiple  solutions,  then  there  are  several 
predictions  for  the  change  of  the  climate;  the  solutions  then  become 
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intransitive.  It  has  been  possible,  for  example,  to  devise  a  rather 
crude  model  of  the  climate  of  the  earth's  surface  that  is  intransi¬ 
tive;  the  model  is  given  by  a  nonlinear,  second  order,  ordinary  dif¬ 
ferential  equation  for  the  earth's  surface  temperature  as  a  function 
of  latitude.  Integrations  of  this  equation  using  the  same  boundary 
conditions  yield  five  climates;  two  of  these  are  physically  unidenti¬ 
fiable,  vrfiile  the  other  three  correspond  to  the  current  climate,  a 
climate  of  the  ice  age,  and  a  climate  for  an  ice  covered  earth,  re¬ 
spectively  (Fig.  8,  Section  IV).  Since  the  two  physically  unidenti¬ 
fiable  solutions  can  be  made  to  vanish  by  the  formulation  of  a  coef¬ 
ficient  of  the  differential  equation,  it  becomes  clear  that  the  use 
of  an  accurate  model  is  a  necessary,  but  not  sufficient,  condition 
for  the  determination  of  the  transitive  or  intransitive  nature  of  a 
prediction  of  any  advertent  change  of  the  climate.  In  the  absence  of 
ergodic  criteria  to  eliminate  all  but  one  of  any  multiple  solutions 
from  a  given  model,  the  answer  to  the  second  subquestion  must  be  de¬ 
termined  from  numerical  investigations  using  an  accurate  model. 

The  third  subquestion,  concerning  the  sensitivity  of  climate  to 
advertent  changes  in  its  key  processes,  can  be  answered  from  numeri¬ 
cal  investigations  using  an  accurate,  transitive  model  for  climate. 

The  fourth  subquestion,  concerning  the  capability  of  current 
technology  to  produce  changes  in  the  climate,  can  be  illustrated  by 
a  potential  mechanism  for  advertent  climate  modification.  Although 
the  sdlar  inclination  is  low  at  high  latitudes,  the  length  of  the 
solar  illumination  during  summer  at  these  latitudes  is  such  that  the 
magnitude  of  the  incoming  solar  radiation  can  be  greater  than  that  at 
low  latitudes  for  a  significant  fraction  of  each  year.  However,  be¬ 
cause  of  the  high  surface  albedo  at  the  high  latitudes,  the  absorbed 
solar  radiation  becomes  significantly  lower  than  that  at  tihe  low  lat¬ 
itudes  (Fig.  9,  Section  VI).  Thus,  a  decrease  in  the  surface  albedo 
accomplished  by  a  large-scale  removal  of  the  snow  and  permanent  ice 
over  the  Arctic  Ocean  could  have  a  significant  impact  on  the  current 
intensity  and  extent  of  the  meridional  cells  at  the  middle  and  high 
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latitudes.  The  use  of  a  transient  thermodynamic  model  for  sea  ice  is 
successful  in  reproducing  theoretically  both  the  thickness  of  the  cur¬ 
rent  ice  layer  over  the  Arctic  Ocean  and  the  temperature  at  its  sur¬ 
face  (Fig.  10,  Section  VI).  Parametric  results  from  this  model  indi¬ 
cate  that  the  use  of  a  suitable  material  dark  enough,  to  decrease  the 
surface  albedo  from  its  current  value  of  about  0.5  to  only  0.4  would 
remove  the  permanent  ice  layer  over  the  Arctic  Ocean  within  a  very  few 
years . 

The  integrations  of  the  foregoing  model  based  on  variable  wave 
and  amplitude  parameters  of  the  baroclinic  disturbances  can  be  applied 
to  the  conditions  of  an  Arctic  Ocean  without  sea  ice  (which  is  also  an 
environment  corresponding  to  the  interglacial  climate  of  the  current 
Pleistocene  Epoch).  This  procedure  would  provide  a  test  of  the  itera¬ 
tive  process  on  the  temperature  at  the  500-mb  level,  which  is  required 
for  the  solution  of  the  problem  of  advertent  climate  modification. 

This  model,  then,  could  be  used  to  seek  answers  for  the  second  and 
third  subquestions  concerning  the  feasibility  of  advertent  modifica¬ 
tion  of  the  climate  of  the  earth’s  surface  at  middle  and  high  northern 
latitudes . 

The  modest  computer  capability  that  is  required  for  the  use  of 
the  model  for  climate ,  based  on  variable  wave  and  amplitude  parameters 
of  *:he  baroclinic  disturbances,  suggests  that  a  world  superpower  can 
solve  the  feasibility  problem  of  advertent  climate  modification  with¬ 
out  the  need  of  new  electronic  computers  of  larger  capabilities. 
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I.  INTRODUCTION 


The  main  objective  of  this  paper  is  to  consider  the  feasibility 
of  controlled  modifications  of  the  climate  of  the  earth’s  surface  at 
middle  and  high  northern  latitudes.  The  concept  of  controlled  or  ad¬ 
vertent  climate  modification  is  different  from  both  the  concept  of 
inadvertent  modifications  of  the  earth’s  climate  by  modern  man’s  ac¬ 
tivities  and  ti>e  concept  of  advertent  weather  modification.  The  con¬ 
cept  of  advertent  climate  modification  stems  from  the  trend  of  world 
superpowers  to  make  timely  exploitations  of  breakthroughs  in  tech¬ 
nology  towards  the  achievement  of  their  economical,  political,  and 
military  goals.  The  considerations  of  the  feasibility  of  advertent 
climate  modification  must  strive  to  gain  a  maximum  understanding  of 
the  key  processes  that  control  atmospheric  and  oceanic  phenomena  near 
the  earth’s  surface. 

The  general  study  of  atmospheric  and  oceanic  processes  during  re¬ 
cent  years  emphasizes  the  use  of  theoretical  models  of  increasing  com¬ 
plexity  that  require  ever  larger  electronic  computers.  This  approach 
yields  a  wealth  of  numerical  details,  which  do  not  necessarily  enhance 
the  understanding  of  the  relevant  controlling  processes.  The  overall 
considerations  of  advertent  climate  modification  have  therefore  at¬ 
tempted  to  capitalize  on  the  availability  of  such  complex  models 
(e.g.,  Rapp,  1970);  however,  as  a  complementary  activity,  a  limited 
manpower  effort  has  also  been  dedicated  to  the  task  of  attempting  to 
gain  insight  into  the  key  processes  that  control  the  climate  at  the 
earth's  surface  through  the  development  of  simpler  models  that  stress 
physical  phenomenology.  Such  simpler  models  provide  two  additional 
advantages  at  a  sacrifice  in  the  wealth  of  numerical  details:  (1) 
they  increase  the  range  and  number  of  parameters  that  must  be  considered 


1 


in  preliminary  numerical  investigations,  a  task  that  would  involve 
prohibitive  investments  of  resources  through  the  use  of  the  more  com¬ 
plex  models;  and  (2)  they  aid  in  the  evaluation  of  the  capability  of 
world  superpowers  that  may  lag  in  the  development  of  ever  larger  elec¬ 
tronic  computers  for  their  consideration  of  the  feasibility  of  adver¬ 
tent  climate  modification.  This  paper,  which  is  based  on  limited  man¬ 
power,  uses  therefore  this  complementary  approach  for  the  considera¬ 
tions  of  the  feasibility  of  advertent  modifications  of  the  climate  of 
the  earth’s  surface  at  middle  and  high  latitudes. 

The  basic  aim  of  emphasizing  the  physical  phenomenology  relevant 
to  the  advertent  modification  of  the  climate  at  the  earth's  surface 
has  shaped  the  scope  of  this  paper,  which  attempts  to  present  the  rel¬ 
evant  physical  and  mathematical  developments  in  sufficient  detail  to 
bring  out  as  clearly  as  possible  the  physical  arguments  and  assump¬ 
tions.  For  this  reason,  the  scope  of  this  paper  includes  the  follow¬ 
ing  topics: 

1.  A  brief  description  of  the  basic  processes  that  control  the 
climate  of  the  earth's  surface. 

2.  The  predictability  of  changes  of  the  climate  at  the  earth's 
surface  or  of  the  long-term  (i.e.,  monthly,  seasonal,  annual) 
statistics  that  result  from  large-scale  wave  phenomena  in 
the  troposphere,  which  control  the  long-term  general  circu¬ 
lation  or  structure  of  the  air  flow  within  the  troposphere. 

3.  The  characteristics  of  the  climate  at  the  earth's  surface 
corresponding  to  the  two  extremes  of  glacial  and  interglacial 
conditions  of  the  current  Pleistocene  Epoch  as  described  by 
their  difference  in  both  the  structure  of  the  general  circu¬ 
lation  and  the  earth's  surface  temperature  as  a  function  of 
latitude. 

4.  A  sample  of  a  potential  mechanism  for  advertent  climate  modi¬ 
fication, 

5.  Description  of  the  general  circulation  for  climate,  including 
a  brief  review  of  the  key  meteorological  developments  rele¬ 
vant  to  the  understanding  of  climate  ever  the  last  30  years. 
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6.  Evaluation  of  available  theoretical  models  of  the  general 
circulation  in  the  troposphere  applicable  to  climate,  by 
using  the  most  recent  reduction  of  experimental  data. 

7.  Generalization  of  the  formulation  of  the  models  so  as  to 
close  the  gap  between  the  theoretical  results  and  the  meas¬ 
urement  of  the  statistics  that  control  the  current  general 
circulation  and  climate  near  the  earth’s  surface. 

The  relevance  of  the  forego!  ig  topics  becomes  clear  from  the  nature  of 
the  key  questions  in  advertent  climate  modification,  which  are  de¬ 
scribed  in  the  next  section. 
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II.  KEY  QUESTIONS  IN  ADVERTENT  CLIMATE  MODIFICATION 

The  fundamental  question  cf  advertent  climate  modification  has  to 
do  with  the  theoretical  evaluation  of  the  feasibility  of  controlled 
modification  of  the  earth's  climate,  in  the  early  1980s,  by  a  world 
superpower  for  either  political  or  military  purposes.  The  1980s  time 
scale  is  fixed  by  the  time  lag  between  initiation  and  operation  of  a 
major  new  project,  which  is  usually  measurable  in  a  large  fraction  of 
a  decade;  and  the  fact  that  it  must  use  proven  and,  therefore,  cur¬ 
rent  advanced  technology.  The  timeliness  of  this  feasibility  ques¬ 
tion  has  been  brought  about  by  two  breakthroughs  in  modem  meteorology: 

1.  Advent  of  meteorological  satellites,  which  provides  a  capa¬ 
bility  for  near  real-time,  worldwide  measurements  of  the 
temperature  distribution  throughout  the  troposphere  (e.g., 
Hidalgo,  1969).  This  capability  includes  measurements  of  the 
wind  field  from  selective  cloud  motion,  and  of  course,  the 
derived  geostrophic  wind  at  middle  and  high  latitudes  from 
the  temperature  field. 

2.  Development  of  modern  electronic  computers,  which  yields  a 
capability  for  both  the  reduction  of  a  large  amount  of  satel¬ 
lite  data  and  the  execution  of  more  refined  theoretical  cal¬ 
culations. 

These  breakthroughs  can  therefore  accelerate  the  progress  towards  a 
more  refined  understanding  of  the  processes  that  sustain  the  general 
circulation  of  the  air  in  the  troposphere,  and  the  subsequent  identi¬ 
fication  of  the  key  processes  for  advertent  climate  modification,  if 
they  exist. 
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The  meteorological  merits  of  the  subject  of  advertent  climate 
modification  are  perhaps  best  addressed  by  quoting  the  late  eminent 
meteorologist,  C.G.  Rossby  (The  Atmosphere  and  Sea  in  Motion,  1959): 

Our  understanding  of  the  general  circulation 
of  the  atmosphere  is  increasing  rapidly,  and  the 
expenses  that  would  go  with  a  large  scale  experi¬ 
ment  of  changing  the  albedo  of  the  earth  do  not 
seem  unsurmountable .  It  is  perhaps  not  surprising 
then  that  some  very  outstanding  scientists  con¬ 
sider  the  question  of  climate  control  to  be  a  real 
and  practical  task  for  meteorological  research. 

As  an  example,  the  highly  judicious  mathematician 
J.  von  Neumann  made  this  statement  in  1955: 
probably  intervention  in  atmospheric  and  climatic 
matters  will  come  in  a  few  decades  and  will  un¬ 
fold  on  a  scale  difficult  to  imagine  at  present. 

Rossby  added: 

Most  meteorologists  would  be  reluctant  perhaps  to 
accept  von  Neumann’s  forecast  today,  but  one 
should  remember  that  during  the  last  decades  the 
technological  development  has  time  after  time 
shown  the  dreams  of  a  visionary  mind  to  be  closer 
to  reality  than  the  commonsense  judgment  of  the 
realists. 

It  may  only  be  appropriate,  more  than  a  decade  later,  to  emphasize 
Rossby’ s  remarks  concerning  the  understanding  of  the  processes  that 
control  the  general  circulation  of  the  atmosphere,  for  the  execution 
of  any  large-scale  experiment  for  the  purpose  of  climate  control  would 
require  an  evaluation  of  the  probable  consequences  to  man’s  environ¬ 
ment  from  the  nonlinear  response  of  the  troposphere  to  such  proposed 
experiment. 

The  fundamental  question  concerning  the  feasibility  of  advertent 
climate  modification  can  be  broken  down  into  the  following  four  basic 
subquestions: 

1.  What  are  the  key  processes  that  control  the  climate  at  the 
earth’s  surface? 

2.  Can  the  results  of  these  key  processes  be  predicted  uniquely 
from  the  use  of  the  physical  laws  controlling  such  processes? 
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3. 


What  is  the  theoretical  sensitivity  of  the  climate  at  the 
earth's  surface  to  changes  in  these  key  processes? 

4.  Does  man  have  the  technology  now  to  use  such  sensitivities 
for  the  control  of  his  climate? 

A  positive  answer  to  the  fundamental  question  concerning  the  feasi¬ 
bility  of  advertent  climate  modification  requires  positive  answers  to 
every  one  of  these  subquestions;  i.e.,  key  processes  must  be  identified 
and  described  through  a  model  based  on  physical  laws,  physical  laws 
must  give  unique  predictions  of  climate  changes,  the  climate  must  have 
an  appropriate  degree  of  sensitivity  to  changes  in  the  key  processes, 
and  man  must  possess  the  technology  necessary  for  the  use  of  such  sensi¬ 
tivities.  While  these  subquestions  are  straightforward  to  pose,  the 
answer  to  the  second  one  is  difficult  because  it  must  invoke  as  yet 
underdeveloped  ergodic  theory  that  has  to  do  with  the  long-term  sta¬ 
tistical  properties  of  solutions  of  differential  equations.  Never¬ 
theless,  this  subquestion  might  be  answerable  with  the  use  of  an  ac¬ 
curate  mathematical  model  of  the  phenomena  that  control  the  current 
long-term  general  circulation,  if  such  model  is  found  to  be  success¬ 
ful  in  describing  uniquely  the  coupling  between  the  main  character¬ 
istics  of  each  climate  during  the  time  the  earth  acquired  its  current 
geographical  distribution  of  oceans  and  continents. 
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III. 


PROCESSES  THAT  CONTROL  THE  EARTH’S  CLIMATE 


The  climate  at  the  earth's  surface  is  characterized  by  the 
monthly,  seasonal,  and  annual  average  magnitude  of  two  thermodynamic 
variables:  the  surface  temperature  (T*)  and  the  humidity  (q)  of  the 
air.  These  variables  are,  however,  coupled  through  the  conservation 
equations  to  the  dynamic  variables  or  two  components  of  the  wind, 
i.e.,  the  zonal  wind  (u)  in  the  longitudinal  direction  (X)  and  the 
meridional  wind  (v)  in  the  latitudinal  direction  (0);  the  vertical 
wind  (w)  in  the  altitude  direction  (z)  is  small,  a  condition  that 
yields  hydrostatic  equilibrium  or  an  exponential  decrease  in  pressure 
(p)  with  increasing  altitude.  The  phenomena  controlling  the  climate 
at  the  earth's  surface  take  place  mainly  in  the  troposphere  or  lowest 
layer  of  the  atmosphere,  which  extends  to  only  about  16  km  from  the 
earth's  surface  at  the  equator  and  to  about  8  km  at  the  poles.  The 
troposphere  is  characterized  by  a  decrease  in  air  temperature  with 
increasing  altitude  (i.e.,  dT/dz  <  0  or  dT/dp  >  0);  this  lowest  layer 
contains  most  of  the  air  in  the  atmosphere,  since  the  air  density  at 
the  tropical  tropopause  (or  upper  limit  of  the  troposphere)  drops  to 
about  15  percent  of  the  value  at  the  earth's  surface.*  Thus,  there 
are  two  immediate  implications  for  climate: 

1.  Only  time  averages  of  the  thermodynamic  and  dynamic  varia¬ 
bles  become  of  interest;  this  is  in  contrast  with  weather 
forecasting,  where  emphasis  must  be  placed  on  the  instan¬ 
taneous  values  of  such  variables. 


* 

The  middle  level  of  the  next  layer,  the  stratosphere  (where 
dT/dz  *  0),  is  of  some  interest  because  the  ozone  concentration 
reaches  a  maximum  there. 
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2.  The  upper  rarified  layers  of  the  atmosphere  (i.e.,  upper 
stratosphere,  mesosphere,  and  thermosphere)  with  their  com¬ 
plicated  chemical  kinetics  and  physics,  are  relatively  unim¬ 
portant  for  the  climate  at  the  earth’s  surface. 

These  two  conditions  allow  the  use  of  significant  mathematical  simpli¬ 
fications  in  the  formulation  of  the  theory  of  the  general  circulation 
for  climate. 

Another  condition  that  allows  a  further  mathematical  simplifica¬ 
tion  of  the  theory  of  the  general  circulation  for  climate  is  the 
rather  remarkable  fact  that  species  present  in  the  troposphere  in 
small  amounts  play  a  dominant  role  in  climate  at  the  earth’s  surface. 
The  most  important  of  these  species  is  water;  near  sea  level  it  is 
present  in  the  air  by  less  than  4  percent  by  volume  in  the  vapor 
phase,  and  less  than  1  percent  in  the  liquid  and  solid  phases  (e.g.. 
Miller  and  Thompson,  1970).  Water  in  its  vapor  and  liqui,  phases  in 
the  lower  levels  of  the  troposphere  influences  the  climate  at  the 
earth’s  surface  mainly  through  its  spectral  properties  for  the  ab¬ 
sorption  of  long-wave  terrestrial  radiation.  Its  spectral  properties 
allow  it  to  absorb  a  fraction  of  the  long-wave  radiation  emitted  by 
the  earth’s  surface,  and  emit  a  downward  flux  of  long-wave  radiation 
back  to  the  earth’s  surface.  This  greenhouse  effect  maintains  the 
temperature  of  the  earth’s  surface  at  a  higher  value  (i.e.,  by  nearly 
50°F,  Willett  and  Sanders,  1959)  than  it  would  be  otherwise.  Water 
also  has  a  significant  impact  on  the  net  heating  of  the  troposphere 
through  its  vaporization-condensation  cycle,  i.e.,  by  absorbing  its 
latent  heat  during  vaporization  at  the  earth’s  surface,  and  releasing 
it  to  heat  the  troposphere  during  its  condensation  vifrtile  ascending  to 
the  middle  and  upper  levels  of  the  troposphere .  The  cloud  formation 
that  results  from  this  condensation  process  has  a  further  effect  on 
the  net  heating  of  the  troposphere,  since  the  statistics  of  the  cloud 
distribution  within  the  troposphere  have  an  important  effect  on  the 
reflection  of  incoming  short-wave  solar  radiation  or  albedo  of  the 
troposphere.  Since  the  vaporization-condensation  cycle  introduces  a 
variable  content  of  water  within  the  troposphere,  the  magnitude  of  the 
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air’s  specific  humidity  as  a  function  of  latitude,  longitude,  and  al¬ 
titude  (or  pressure,  since  these  two  variables  are  interchangeable 
from  the  condition  of  hydrostatic  equilibrium  in  the  vertical  direc¬ 
tion)  must  be  determined  from  the  conservation  equation  for  water 
vapor.  However,  the  small  concentration  of  water  vapor  relative  to 
the  air  allows  a  mathematical  decoupling  of  the  equations  of  motion 
from  the  continuity  equation  for  the  specific  humidity.* 

Finally,  the  rotation  of  the  earth  and  the  orientation  of  its 
axis  of  rotation  relative  to  the  sun  allow  the  use  of  zonal  (or  longi¬ 
tudinal)  averages  for  the  thermodynamic  and  dynamic  variables  of  in¬ 
terest  for  climate;  i.e.,  the  rotation  of  the  earth  tends  to  average 
the  magnitude  of  such  variables  at  a  given  latitude  whenever  the  time 
interval  of  interest  becomes  appreciably  longer  than  one  day.  The 
zonal  averaging  allows  the  value  of  a  dependent  variable  at  any  longi¬ 
tude  to  be  expressed  as  the  stun  of  its  zonal  average  and  the  deviation 
from  its  zonal  average.  Tho  former  is  usually  referred  to  as  pertain¬ 
ing  to  the  "circulation,”  and  the  latter  to  the  "eddy,"  or  large-scale, 
two-dimensional  turbulent  motion.  Since  longitude  and  time  are  inde¬ 
pendent  variables  of  the  atmospheric  motion,  a  time-average  of  a  de¬ 
pendent  variable  at  a  given  longitude  also  yields  two  components,  i.e., 
the  long-term  time  average  (monthly,  seasonal,  or  annual)  or  "stand¬ 
ing"  motion  and  its  deviation  from  the  long-term  time  average  or 
"transient"  motion.  A  dependent  variable  at  any  longitude  and  time 
is  then  expressible  by  the  sum  of  four  components:  the  standing  and 
transient  circulations  and  the  standing  and  transient  eddy  motions. 
Thus,  the  meridional  wind  v,  for  example,  is  given  by 

v  =  [v]  +  [v]#  +  v*  +  v*  ' 


Another  specie  of  interest  is  carbon  dioxide.  However,  since  carbon 
dioxide  is  well  mixed  with  air,  its  mixing  ratio  (3  x  10“4)  is  con¬ 
stant. 
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where  the  bracket,  star,  bar,  and  prime  denote,  respectively,  the 
zonal  average,  eddy,  time  average,  and  transient  components.  Typical 
examples  of  these  types  of  motion  are  as  follows  (Lorenz,  1967): 

1.  Trade  winds  at  low  latitudes  for  the  standing  circulation. 

2.  Asiatic  summer  and  winter  monsoons,  or  features  in  addition 
to  those  of  the  standing  circulation  which  appeal  vdien  the 
variables  are  averaged .with  respect  to  time  alone. 

3.  Fluctuations  of  the  zonal  index,  or  features,  in  addition  to 
those  of  the  standing  circulation  which  appear  when  the  vari¬ 
ables  are  averaged  with  respect  to  longitude  alone. 

4.  Migratory  cyclones,  or  features,  in  addition  to  those  of  the 
standing  circulation  which  appear  when  the  variables  are  not 
averaged. 

The  long-term  zonal  average  of  the  poleward  transport  of  any  quantity 
X  (where  X  denotes  any  quantity  such  as  a  velocity  component,  tempera¬ 
ture,  specific  humidity,  etc.)  may  then  be  resolved  into  the  amounts 
accomplished  by  the  separate  components  of  v,  i.e., 

[Xv  ]  =  [X][v]  +  [xrtv]'  +  [X*v*  ]  +  [X*  'v* '  ] 


The  terms  at  the  right-hand  side  then  represent,  respectively,  the 
standing  circulation  or  cell  transport,  transient  cell  transport, 
standing  eddy  transport,  and  the  transient  eddy  transport. 

The  foregoing  considerations  emphasize  general  characteristics 
of  atmospheric  and  planetary  phenomena  that  are  relevant  to  the  formu¬ 
lation  of  a  theory  of  the  general  circulation  for  climate  applica¬ 
tions.  The  specific  description  of  the  key  processes  that  control 
the  climate  at  the  earth’s  surface  involves  two  distinct  but  coupled 
types  of  phenomena: 

1.  The  net  heating  rate  of  the  atmosphere,  vdiich  is  a  result  of 
the  following  four  processes: 
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a.  Absorption  of  long-wave  terrestrial  radiation  mainly  by 
water  and  carbon  dioxide  in  the  troposphere. 

b.  Absorption  of  short-wave  solar  radiation  mainly  by  oxygen 
and  ozone  in  the  stratosphere;  the  absorption  of  solar 
radiation  by  these  species  is  relatively  small,  and  the 
atmosphere  as  a  whole  may  be  assumed  to  be  transparent 

to  the  short-wave  solar  radiation  (Charney,  1959). 

c.  Condensation  of  water  vapor  mainly  in  the  troposphere. 

d.  Flux  of  sensible  heat  from  the  earth’s  surface  trans¬ 
ported  by  small-scale,  three-dimensional  turbulence  in 
the  surface  boundary  layer  or  layer  extending  from  the 
earth’s  surface  to  an  altitude  of  about  1  km. 

2.  Large-scale  poleward  transports  of  angular  momentum,  sen¬ 
sible  heat,  and  water  vapor  by  the  transient  and  standing 
eddy  motion  at  the  middle  and  upper  levels  of  the  tropo¬ 
sphere;  this  motion  is  triggered  by  the  net  heating  rate  of 
the  troposphere.  The  large-scale  transports  of  angular  mo¬ 
mentum  and  sensible  heat  by  the  eddy  motion  then  set  up  the 
dynamics  of  the  general  circulation  in  distinct  cells  of 
meridional  circulation  of  a  statistical  nature. 

Each  of  these  two  types  of  phenomena  is  described  below  in  more  de¬ 
tail. 

A.  NET  HEATING  RATE  OF  THE  ATMOSPHERE 

The  absorption  of  long -wave  terrestrial  radiation  and  short-wave 
solar  radiation  by  water,  carbon  dioxide,  oxygen,  and  ozone  is  con¬ 
trolled  by  the  absorptivity  of  these  species  as  a  function  of  wave 
length  (their  spectra  are  shown,  for  example,  in  Miller-Thompson, 

1970).  The  condensation  of  water  vapor  is  coupled  with  the  dynamics 
of  the  circulation  and  eddy  motions  in  the  troposphere  ( Smagorinsky , 
1963),  while  the  transport  of  sensible  heat  in  the  surface  boundary 
layer  is  governed  by  the  appropriate  energy  equation.  The  instability 
of  the  system  is  controlled  by  the  Richardson  number,  i.e.,  a  criterion 
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that  expresses  the  relative  role  of  the  buoyancy  e’nd  wind  shear  in 
three-dimensional  turbulence  (e.g.,  Hidalgo,  1970;  Monin,  1970).  A 
parameterization  of  these  complex  phenomena  has  been  given  by 
Smagorinsky  (1963),  using  previous  theoretical  and  experimental  re¬ 
sults  by  Phillips  (1956),  Charney  (1959),  Houghton  (1954),  Budyko 
(1956),  and  London  (1957). 

This  parameterization  is  based  on  long-term  (annual)  zonal  aver¬ 
ages  of  the  relevant  parameters  and  is  described  schematically  in  Fig. 
1,  which  shows  three  altitude  or  pressure  levels  at  any  given  lati¬ 
tude:  (1)  top  of  the  atmosphere;  (2)  the  atmosphere,  as  represented 
by  a  column  with  a  mean  temperature  corresponding  to  the  500-mb  pres¬ 
sure  level  or  altitude  of  about  5.6  km  lor  a  standard  atmosphere;  and 
(3)  the  earth’s  surface. 


SHORT-WAVE  LONG-WAVE  NON  RADIATIVE 

SOLAR  RADIATION  TERRESTRIAL  RADIATION  PROCESSES 


FIGURE  1.  Schematic  Diagram  of  the  Processes  Contributing  to  the  Heat  Balance 
in  the  Troposphere  (Smagorinsky,  1963) 
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Figure  1  is  self-explanatory;  at  each  level,  open  circles  denote 
either  emission  or  flux  divergence  and  closed  circles  signify  either 
absorption  or  flux  convergence.  The  left  side  of  the  figure  describes 
the  absorption  of  short-wave  radiation  by  the  atmosphere*  (S  )  and 
earth’s  surface  (S*),  where  SQ  denotes  the  solar  radiation  at  the  top 
of  the  atmosphere,  A  the  albedo  (or  reflection  to  space  of  short-wave 
radiation)  of  the  atmosphere,  X  the  short-wave  absorptivity  by  oxygen 
and  ozone  in  tne  atmospheric  column,  A*  the  albedo  of  the  earth’s  sur¬ 
face,  and  Ap  the  planetary  albedo.  The  middle  of  the  figure  describes 
the  emission  and  absorption  of  long-wave  terrestrial  radiation,  where 
T,v  denotes  the  temperature  at  the  earth's  surface,  o  the  Stephan- 
Boltzmann  constant,  T  the  long-wave  absorptivity  by  the  water  and  car¬ 
bon  dioxide  in  the  atmospheric  column,  and  the  empirical  parameters 
vt  and  v l  pertain  to  the  upward  and  downward  fluxes  of  radiation,  re¬ 
spectively.  The  figure  indicates  that  the  earth’s  surface  absorbs  and 
emits  as  a  blackbody,  while  the  atmosphere  absorbs  and  emits  as  a  gray 
body.  Finally,  the  right-hand  side  of  the  figure  shows  the  nonradia- 
tive  processes. 

As  shown  in  Fig.  1,  for  the  earth’s  surface  level,  flux  diver¬ 
gence  of  heat  in  the  oceans  due  to  lateral  ocean  transports  is  indi¬ 
cated  by  Mq,  and  the  flux  divergence  of  both  sensible  and  latent  heat 
due  to  the  eddy  transports  of  temperature  and  humidity  from  the 
earth’s  surface  to  the  troposphere  is  indicated  by  Eg  and  E^>  respec¬ 
tively.  At  the  500-mb  level,  the  heating  of  the  atmospheric  column 
is  indicated  by  both  the  condensation  of  water  vapor  (C)  and  the  heat 
flux  Eg  from  the  surface  boundary  layer. 

The  net  average  heating  rate  of  the  atmospheric  column  may  be  de¬ 
fine  j  as  Qj,  i.e. , 


*2 


The  notation  [”]  has  been  dropped  to  simplify  this  figure. 


15 


where  Q  is  the  net  heating  rate  at  any  pressure  level,  ami  the  limits 
of  integration  correspond  to  the  top  of  the  atmosphere  (p  =  0)  and  to 
the  earth’s  surface  (p  =  =  1000  mb).  With  this  definition  and  the 

use  of  the  principle  of  conservation  of  energy,  the  closed  and  open 
circles  at  500  mb  of  Fig.  1  yield  the  net  heating  rate  of  the  atmos¬ 
pheric  column  as  follows: 


P4H0  444 

=  Sa  +  roT*  -  vtoTj  -  virfTj  +  c  +  Eg  (1) 

where  the  terms  in  the  right-hand  side  are  obtained  from  annual  ex¬ 
perimental  data. 

Figure  2  shows  the  magnitude  of  each  of  the  terns  in  Eq.  1  for 
the  current  climate,  and  the  rather  dominant  role  of  the  long-wave 
radiation  terms.  If  it  is  assumed  that  there  is  no  heat  accumulation 
at  the  earth’s  surface,  the  circles  corresponding  to  the  earth’s  sur¬ 
face  yield 


ot£  =  S*  +  vidTj  -  (Eg  +  El  +  MQ)  (2) 

i.e.,  the  surface  temperature  is  a  function  of  the  short-wave  radia¬ 
tion  reaching  the  surface,  its  albedo,  the  downward  flux  of  long -wave 
radiation,  and  the  divergence  of  the  sensible,  latent,  and  oceanic 
heat  fluxes.  For  the  conditions  indicated  in  Fig.  2  for  the  current 
climate,  the  surface  and  mean  atmospheric  temperatures  vary  as  a  func¬ 
tion  of  latitude  in  the  range  260  ST*  s  300°K  and  240  s  Tg  £  270°K. 
Figure  3  shows  the  magnitude  of  each  of  the  terms  in  Eq.  2  for  the 
current  climate,  and  it  is  consistent  with  the  annual  mean  data  of  the 
previous  figure,  i.e.,  Budyko's  experimental  values  for  the  lateral 
oceanic  transport  MQ  have  been  reduced  by  9.4  to  make  them  consistent 
with  London’s  data  (Smagorinsky ,  1963).  Figure  3  shows  the  dominant 
effect  of  the  long-wave  radiation  term  in  the  right-hand  side  of  Eq.  2 
at  all  latitudes,  and  the  relatively  small  magnitude  of  the  oceanic 
transport  MQ. 
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NORTH  LATITUDE,  d«9 

FIGURE  2.  Net  Heating  Rate  Per  Unit  Area  of  a  Column,  p^Q^/g,  Based  on 
Annual  Mean  Data  (Smagorinsky,  1963) 
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NORTH  LATITUDE,  dag 


FIGURE  3.  Contribution  to  Earth'*  Surface  Temperature  (Tj  of  Long-Wave  Downward 
Radiation  (Li),  Short-Wave  Radiation  (Sj,  Flux  Divergence  of  Heat  Due 
to  Eddy  Transports  of  Sensible  and  Latent  Heat  (E^  +  E^)  and  Lateral 
Oceanic  Transport  (M  ) 


B.  LARGE-SCALE  POLEWARD  TRANSPORTS 

The  turbulent  nature  of  the  eddy  motion  triggered  by  the  net 
heading  rates  of  the  troposphere  is  demonstrated  by  the  longitudinal 
spectrum  of  kinetic  energy  at  synoptic  wavelengths.  Figure  4  shows 
the  spectral  distribution  of  large-scale  atmospheric  kinetic  energy 
as  a  function  of  pressure  level  and  latitude  (Julian  et  al.,  1970). 

The  kinetic  energy  spectrum  is  defined  as  the  portion  of  the  spectrum 
arising  from  stochastic  transient  motion  on  all  scales,  and  it  is  de¬ 
rived  from  both  reduced  wind  data  and  untreated  observed  winds.  The 
hemispheric  wave  number  (k)  denotes  the  number  of  wavelengths  (X)  in 
a  latitudinal  circle  with  radius  a*cos  0,  where  a  is  the  radius  of 
the  earth  and  0  the  latitude  of  the  circle,  i.e., 

,  2na  cos  0 

k  X 

Thus,  large  wavelengths  or  large-scale  phenomena  correspond  to  small 
wave  numbers.  The  kinetic  energy  in  Fig.  4  has  been  normalized  with 
respect  to  its  value  at  k  =  6.  The  two-dimensional  character  of  the 
turbulent  motion  in  the  troposphere  is  indicated  in  Fig.  4  by  the 
spectral  slope  at  the  higher  wave  numbers;  i.e.,  all  estimates  of  the 
spectrum  suggest  that  the  spectral  slope  in  the  3500-  to  1500-km  wave¬ 
length  region  can  be  represented  by  a  power  law  with  an  exponent  in 
the  range  -2.7  to  -3.0,  which  is  different  tnan  the  -5/3  value  for 
three-dimensional  or  isotropic  turbulence  at  very  small  wavelengths. 
The  two-dimensional  nature  of  the  turbulence  at  the  high  wave  numbers 
becomes  important  for  considerations  of  the  theoretical  limit  on  long- 
range  weather  forecasting,  vrtiich  is  estimated  to  be  at  most  no  more 
than  a  few  weeks.  For  climate,  however,  the  interest  is  on  large- 
scale  phenomena  or  low  wave  numbers,  since  a  theoretical  solution  of 
barotropic  atmospheric  motion  (involving  the  use  of  the  momentum  and 
continuity  equations)  and  a  homogeneous  earth  indicate  that  the  wave 
number  for  a  large-scale  stationary  disturbance  is  given  by  (Haurwitz, 
1940) 
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where  [u]  denotes  the  standing  zonal  circulation.  For  example,  this 
barotropic  solution  yields  k  %  6  when  [u]  =  20  m/sec  at  45-deg  lati¬ 
tude.  As  shown  later,  the  effect  of  baroclinic  disturbances  (where 
the  energy  equation  is  used  together  with  the  momentum  and  continuity 
equations)  and  a  nonhomogeneous  earth's  surface  (as  given  by  the  dis¬ 
tribution  of  oceans  and  continents)  may  yield  a  wave  number  appropriate 
for  climate  that  is  even  lower  than  the  above  barotropic  value.  A 
fundamental  difference  of  interest  for  climate  between  large-scale 
two-dimensional  turbulence  at  very  low  wave  numbers,  and  small-scale 
three-dimensional  turbulence  at  very  high  wave  numbers,  is  the  fact 
that  the  transports  by  the  former  can  be  countergradient,  i.e.,  the 
usual  concepts  of  eddy  viscosity  for  very  small  or  microscale  turbu¬ 
lence  cannot  be  extended  to  the  large-scale  transports.  The  eddy  vis¬ 
cosity  or  mixing-length  concepts  in  microturbulence  assume  that  the 
transports  are  given  by  the  product  of  an  eddy  viscosity  and  the  gra¬ 
dient  of  the  average  flow.  This  assumption  breaks  down  for  the  large- 
scale  transport  of  angular  momentum  at  middle  and  high  latitudes,  for 
example,  since  at  these  latitudes  the  transport  is  positive  while  the 
latitudinal  gradient  of  the  average  flow  (i.e.,  d[u]/d6)  is  negative. 
Since  this  transport  and  the  gradient  d[u]/d9  remain  both  positive  at 
lower  latitudes,  even  a  concept  of  a  "negative”  eddy  viscosity  becomes 
unapplicable  at  any  arbitrary  latitude. 

The  organization  of  general  circulation  dynamics  into  distinct 
cells  by  the  large-scale  transports  of  angular  momentum  and  sensible 
heat  is  illustrated  in  Fig.  5  for  the  current  climate.  The  figure 
shows  three  meridional  cells:  (1)  the  direct  Hadley  cell,  vthere  humid 
air  ascends  near  the  equator  and  dry  air  descends  near  the  Horse  lati¬ 
tude;  (2)  the  indirect  Ferrel  cell,  where  dry  air  descends  near  the 
Horse  latitude  and  moist  air  ascends  at  the  higher  latitude;  and  (3) 
the  polar  cell  at  the  high  latitudes.  The  Hadley  cell  is  controlled 
by  the  poleward  large-scale  transport  of  angular  momentum,  which  reaches 
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a  maximum  at  a  latitude  of  about  30  deg  and  near  the  250-mb  pressure 
level  (e.g.,  see  Lorenz,  1967)  or  at  an  altitude  of  about  10  km  for  a 
standard  atmosphere,  an  altitude  that  corresponds  to  that  of  the  jet 
streams.  The  convergence  of  the  eddy  transports  of  angular  momentum 
acts  as  a  mechanical  force,  which  induces  the  downward  motion  of  the 
meridional  cells  at  the  30-deg  latitude.  Similarly,  the  transport  of 
sensible  heat  reaches  its  maximum  value  near  the  50-deg  latitude;  the 
convergence  of  the  eddy  transport  of  sensible  heat  acts  as  a  heat 
source,  which  induces  the  upward  motion  of  the  indirect  Ferrel  and 
polar  cells  near  that  latitude.  The  opposite  directions  of  the  meridi¬ 
onal  winds  in  the  lower  branches  of  the  Hadley  and  Ferrel  cells,  to¬ 
gether  with  the  effect  of  the  Corioles  acceleration  (-25  x  [v])  from 
the  earth’s  rotation  (5),  control  the  direction  of  the  zonally  aver¬ 
aged  standing  winds  at  the  earth’s  surface,  which  are  as  shown  in 
Fig.  5. 

Figure  6  shows  the  annual  average  precipitation  and  evaporation 
as  a  function  of  latitude  for  the  current  climate.  Precipitation 
reaches  a  maximum  near  the  equator  or  the  ascending  branch  of  the  di¬ 
rect  Hadley  cell  in  the  northern  hemisphere.  This  maximum  precipita¬ 
tion  yields  a  corresponding  maximum  at  the  same  latitude  in  the  net 
heating  of  the  troposphere  by  the  condensation  process  (i.e.,  as  im¬ 
plied  in  Fig.  2).  Figure  6  also  shows  that  the  precipitation  reaches 
a  second  peak  near  50  deg  or  near  the  latitude  of  the  ascending  branch 
of  the  Ferrel  cell.  Figure  2  again  indicates  a  corresponding  second 
peak  at  the  same  latitude  in  the  net  heating  by  condensation.  The 
maximum  in  precipitation  near  the  equator,  where  the  NE  and  SE  sur¬ 
face  winds  converge  (Fig.  5),  is  referred  to  as  the  intertropical  con¬ 
vergence  zone  ( ITCZ )  or  equatorial  doldrums.  Note  that  Fig.  6  does 
not  show  a  double  ITCZ  across  the  equator,  as  would  be  expected  from 
the  Hadley  cells  in  each  hemisphere  (Fig.  5).  A  numerical  model  based 
on  a  troposphere  interacting  with  the  oceans  has  been  successful  in 
reproducing  the  single  ITCZ  north  of  the  equator  (Pike,  1970). 
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FIGURE  6.  Average  Annual  Evaporation  and  Precipitation  (Lorenz,  1967) 

Figure  ?  shows  results  from  this  model  for  the  precipitation  and 
ocean  surface  temperature  as  a  function  of  time  during  the  88  days 
used  in  these  calculations.  The  dashed  lines  show  the  ITCZ  and  the 
temperature  during  the  first  nine  days,  vtfien  the  ocean’s  surface  tem¬ 
perature  remains  constant;  these  lines  show  a  single  ITCZ,  which  re¬ 
mains  fixed  almost  at  the  equator.  The  westward  surface  winds  induce 
a  westward  sea  current  at  the  equator,  while  the  Coriolis  acceleration 
produces  poleward  (i.e.,  opposite)  meridional  sea  currents  in  each 
hemisphere  at  the  equator.  The  result  is  an  entrainment  of  deeper, 
colder  water  at  the  equator;  in  these  calculations,  this  effect  starts 
to  cool  the  sea  surface  temperature  after  the  ninth  day,  and  it  be¬ 
comes  about  3  deg  colder  after  88  days  (Fig.  7).  As  the  temperature 
starts  to  cool  after  the  ninth  day,  the  ITCZ  moves  northward  slowly; 
after  88  days  it  reaches  6  deg  N  over  a  predicted  surface  temperature 
maximum.  Again,  the  model  does  not  yielc  a  split,  in  the  ITCZ,  into 
two  parts,  each  centered  over  a  temperature  maximum.  Satellite  data 
for  1965  through  1967  indicate  that  the  zonal  average  of  the  latitude 
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The  foregoing  considerations  indicate  that,  for  the  current  cli¬ 
mate,  the  Hadley  cell  controls  tropical  meteorological  phenomena, 
while  the  Ferrel  and  polar  cells  control  these  phenomena  at  the  middle 
and  high  latitudes.  This  implies  that  it  may  be  feasible  to  decouple 
tropical  phenomena  from  considerations  of  the  effect  of  advertent  cli¬ 
mate  modification  at  the  middle  and  high  latitudes  on  the  structure 
of  the  current  Ferrel  and  polar  cells. 

The  description  of  the  key  processes  controlling  the  current  cli¬ 
mate  at  the  earth's  surface  identifies  the  phenomena  of  large-scale, 
two-dimensional  turbulence  at  the  middle  and  upper  levels  of  the  tro¬ 
posphere  as  being  of  fundamental  importance.  Such  phenomena  con¬ 
trol  the  latitudinal  variation  of  the  mean  tropospheric  temperature 
(T2)  through  the  large-scale  transport  of  sensible  heat,  the  current 
configuration  of  the  three  cells  in  the  troposphere,  and  therefore, 
the  distribution  of  specific  humidity.  The  mean  tropospheric  tempera¬ 
ture  and  humidity  control  the  downward  flux  of  long-wave  radiation 
(i.e.,  vloT^,  Eq.  2),  which  has  a  dominant  effect  on  the  other  basic 
variable  of  advertent  climate  modification,  i.e.,  T*(0),  as  shown  in 
Fig.  3. 


IV.  PREDICTABILITY  OF  THE  CLIMATE  AT  THE  EARTH’S  SURFACE 


As  shown  later,  the  parameterization  of  the  net  heating  rates  of 
the  troposphere  (i.e.,  Eq.  1)  can  be  used,  together  with  the  conserva¬ 
tion  equations  and  the  characteristics  of  large-scale  wave  phenomena, 
to  determine  the  statistics  that  control  the  cell  structures  of  the 
current  long-term  general  circulation.  These  statistics  are,  for  ex¬ 
ample,  the  annual,  large-scale  transports  of  angular  momentum  (at  250 
mb),  sensible  heat  and  humidity  (at  500  mb),  as  well  as  the  annual 
meridional  and  zonal  kinetic  energies  and  circulations.  This  mathe¬ 
matical  description  of  the  key  processes  becomes  a  model  that  may  be 
used  for  predicting  changes  in  the  long-term,  zonal  averager,  of  the 
surface  temperature,  T*(0),  and  humidity  from  advertent  (or  inadver¬ 
tent  for  that  matter)  climate  modifications.  The  second  basic  sub¬ 
question  identified  earlier  arises  then  from  considerations  of  the 
uniqueness  (or  lack  of  it)  of  the  predicted  T*(  0)  that  is  derived, 
for  example,  from  such  a  nonlinear  mathematical  model. 

There  are  two  different  aspects  concerning  the  general  question 
of  the  predictability  of  tropospheric  phenomena.  One  is  short¬ 
term,  relevant  to  the  transient  phenomena  of  interest  for  weather 
forecasting  (e.g.,  Thompson,  1957;  Smagorinsky,  1963;  Lorenz,  1968; 
etc.).  The  other  is  long-term,  relevant  to  the  statistical  properties 
of  the  state  of  the  troposphere  of  interest  for  climate  (e.g.,  Lorenz, 
1967,  1968,  1970).  The  former  is  concerned  with  the  effect  of  error 
magnitude  in  the  initial  conditions  on  the  results  of  short-term  nu¬ 
merical  integrations  of  the  conservation  equations.  The  latter  in¬ 
quires  into  the  uniqueness  of  the  long-term  statistics  of  the  solu¬ 
tions.  For  climate,  then,  the  nonlinear  models  describing  the  key 
processes  can  be  either  (1)  transitive,  which  yield  for  a  given  set  of 
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condition?"  only  one  set  of  long-term  statistics  or  one  climate;  and 
(2)  intransitive,  which  yield  again  for  a  given  set  of  conditions  two 
or  more  sets  of  statistics;  consequently,  each  set  of  statistics  be¬ 
comes  an  alternative,  physically  possible  climate.  Since  there  is  no 
criterion  based  on  physical  law  for  the  selection  of  only  one  of  these 
alternative  climates,  it  is  not  currently  known  whether  the  atmosphere- 
ocsan-earth  system  is  transitive  or  intransitive.  In  this  regard,  it 
is  of  theoretical  interest  to  consider  the  following  example  of  an  in¬ 
transitive  nonlinear  equation,  vhich  is  a  crude  approximation  of  the 
climates  description,  i.e.,  T*(0),  at  the  earth's  surface. 

The  first  law  of  thermodynamics  may  be  applied  to  a  control  vol¬ 
ume  (e.g.,  Shapiro,  1953)  for  the  earth-atmo sphere -ocean  system,  i.e., 

2 

a  volume  defined  at  a  latitudinal  circle  by  an  area  element  2na  cosede 
on  the  earth's  surface  and  a  unit  height  of  the  planetary  boundary 
layer.  For  quasi-steady  conditions,  the  net  flux  of  energy  normal  to 
the  earth's  surface  must  equal  the  net  poleward  flow  of  energy  (Ep^v^e^) 
across  the  area  2TTacos0  •  1;  hence,  using  the  notation  of  Fig.  1,  one 
obtains  (Faegre,  1972): 


The  poleward  flow  of  energy  consists  of  three  components:  the  trans¬ 
ports  of  sensible  and  latent  heat  in  the  atmosphere  and  the  transport 
of  sensible  heat  in  the  top  layer  of  the  ocean.  Recalling  that  there 
is  no  net  flux  of  energy  in  a  normal  direction  at  the  earth's  surface 
(i.e.,  Eq.  2),  the  right-hand  side  of  Eq.  3  denotes  the  net  energy 
flux  per  unit  area  at  the  top  of  the  atmosphere.*  By  using  the  crude 

4 

assumptions  that  Lt  =  KoT*,  and  that  each  of  the  three  energy  terms 
in  the  left-hand  side  of  Eq.  3  is  proportional  to  dT*/de,  Eq.  3  becomes 


An  alternative  method  to  derive  Eq.  3  is  the  use  of  a  procedure 
similar  to  that  of  Eq.  6. 
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a  second-order,  ordinary,  nonlinear  differential  equation  for  de¬ 
scribing  the  climate  at  the  earth's  surface  through  T*(  0).  The  two 
required  boundary  conditions  must  specify  that  there  is  no  energy  flow 
off  either  pole,  i.e.,  dT*/d0  -•  0  as  0  approaches  90  deg  N  and  90 
deg  S. 

Figure  8  shows  the  results  of  numerical  integrations  of  Eq.  3 
using  present  values  of  the  parameters  in  this  model;  the  figure  shows 
that  there  are  five  solutions  to  Eq.  3,  or  five  climates.  The  solid 
line  corresponds  to  the  present  climate,  the  higher  dashed  line  to  a 
climate  that  might  correspond  to  that  of  an  ice  age,  the  lower  dashed 
line  to  a  climate  for  an  ice -covered  earth,  and  the  two  dotted  lines 
to  physically  unidentifiable  climates.  The  two  assumptions  used  in 
Eq.  3  have  a  major  impact  on  these  solutions,  i.e.,  (1)  if  Lt  is  for¬ 
mulated  with  a  constant  value  of  K  corresponding  to  the  appropriate 
average  value  with  latitude,  the  two  dotted  solutions  in  Fig.  8  disap¬ 
pear;  and  (2)  if  the  formulation  of  the  transports  are  modified  to  in¬ 
clude  effects  of  the  meridional  circulations,  then  the  model  using  a 
variable  K  becomes  transitive  and  similar  to  other  crude  models  for 
the  climate  at  the  earth's  surface  (i.e.,  Budyko,  1969,  1970;  Sellers, 
1969,  1970).  Thus,  the  transitive  or  intransitive  nature  of  a  model 
becomes  sensitive  to  the  built-in  structure  of  such  a  nonlinear  model. 
Stated  in  another  way,  the  use  of  a  most  accurate  structure  in  a  non¬ 
linear  model  becomes  a  necessary,  but  not  sufficient,  condition  for 
defining  either  its  transitive  or  intransitive  nature.  In  the  absence 
of  ergodic  criteria  to  eliminate  all  but  one  of  the  multiple  solutions 
for  a  given  set  of  conditions,  the  transitive  or  intransitive  nature 
of  a  nonlinear  model  with  a  most  accurate  structure  must  depend  on  the 
results  from  investigations  aimed  at  showing  an  absence  of  more  than 
one  solution. 
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FIGURE  8.  Surface  Temperature  as  a  Function  of  Latitude  Predicted  from  the 
Same  Mathematical  Formulation,  Eq.  3  (Faegre,  1972) 
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V.  CHARACTERISTICS  OF  THE  CLIMATES  OF  THE 
PLEISTOCENE  EPOCH 

The  transitive  or  intransitive  nature  of  nonlinear  models  for  de¬ 
scribing  changes  of  climate  at  the  earth's  surface  leads  to  an  inter¬ 
est  in  the  main  characteristics  of  the  two  extreme  types  of  climate 
which  have  existed  from  the  time  the  earth  acquired  its  current  geo¬ 
graphical  distribution  of  oceans  and  continents,  i.e.,  the  major 
glacial-interglacial  cycles  of  climate  and  atmospheric  circulation  of 
the  current  Pleistocene  Epoch.  The  problem  of  the  changes  from  in¬ 
terglacial  to  glacial  climates  is  different  from  that  of  advertent  or 
inadvertent  climate  modifications.  This  is  because  the  former  has  an 
additional  degree  of  freedom  through  the  variable  magnitude  of  the  in¬ 
coming  solar  radiation,  SQ(0),  as  caused  by  changes  in  both  solar  and 
earth's  orbital  characteristics  over  extremely  long  times.  Neverthe¬ 
less,  the  characteristics  of  the  deduced  general  circulations  corre¬ 
sponding  to  the  glacial  and  interglacial  extremes  may  be  used  to  test 
the  transitive  or  intransitive  nature  of  a  model  that  describes  the 
current  climate  accurately  and  uniquely.  This  test  would  involve  an 
evaluation  of  the  capability  of  the  model  to  reproduce  transitively 
the  coupling  between  the  deduced  cell  structure  of  the  general  circu¬ 
lation  and  the  net  heating  of  the  troposphere  under  glacial  or  inter¬ 
glacial  conditions.  Note  that  such  a  test  would  avoid  the  difficult 
problem  concerning  a  prediction  of  the  net  heating  of  the  troposphere 
under  glacial  conditions,  i.e.,  the  rather  controversial  predictions 
of  the  onset  and  recess  of  ice-age  phenomena  (e.g.,  Mitchell,  1968; 
MacCracken,  1969;  Schell,  1971). 

The  end  of  the  most  recent  (Wtfrm)  major  glacial  period  of  the 
Pleistocene  Epoch  is  usually  placed  at  approximately  6500  B.C.  It 
was  the  time  at  which  the  ice  boundaries  and  climate  conditions  in 
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Europe  had  returned  essentially  to  those  of  today.  The  probable  main 
characteristics  of  the  general  circulation  and  climate  during  the 
Pleistocene  glacial-interglacial  cycles  have  been  deduced  from  geolog¬ 
ical  evidence  (e.g.,  Willett,  1953).  These  characteristics  indicate 
that  the  statistical  cell  structure  of  the  meridional  circulation 
(i.e.,  Fig.  5)  persists  even  under  the  extremes  of  the  minimum  and 
maximum  glaciations  of  the  Pleistocene  Epoch.  The  available  evidence 
indicates  that  these  major  glacial-interglacial  cycles  are  in  phase 
rather  than  in  opposition  in  the  earth's  northern  and  southern  hemi¬ 
spheres.  The  general  circulation  for  the  conditions  of  minimum  and 
maximum  glaciation  is  characterized,  respectively,  by  two  and  three 
meridional  cells,  i.e.,  by  the  configuration  of  the  polar  cell,  which 
disappears  for  minimum  glaciation  but  for  maximum  glaciation  extends 
equatorward  by  about  15  deg  from  the  latitude  of  its  current  ascending 
branch  (e.g.,  Figs.  5,  6).  Thus,  for  minimum  glaciation,  the  surface 
easterlies  at  high  latitudes  disappear  with  the  polar  cell,  while  the 
westerlies  of  the  Ferrel  cell  expand  polewaxd  to  the  poles.  The  pole- 
ward  expansion  of  the  westerlies  induces  a  poleward  expansion  of  the 
Hadley  cell,  i.e.,  the  latitude  of  the  descending  branches  of  the 
Ferrel  and  Hadley  cells  increases  from  its  current  value.  Conversely, 
for  maximum  glaciation,  the  easterlies  of  the  polar  cell  are  displaced 
towards  the  equator;  this  effect  induces  a  displacement  of  the  de¬ 
scending  branch  of  the  Hadley  cell  to  latitudes  lower  than  the  cur¬ 
rent  value.  For  conditions  of  maximum  glaciation,  the  Hadley  cell  is 
characterized  by  an  ITCZ  that  is  stronger  than  the  current  maximum 
precipitation  near  the  equator. 

The  climate  corresponding  to  the  two-  and  three-cell  structures 
of  the  minimum  and  maximum  glaciations  is  as  follows: 

1.  For  interglacial  conditions,  the  ascending  branch  of  the 
Ferrel  cell  at  the  poles  yields  an  increase  of  warmth  and 
precipitation  relative  to  current  values  at  the  polar  lati¬ 
tudes;  the  result  is  an  absence  of  permanent  glaciation  on 
the  earth's  surface  at  these  high  latitudes  and  a  displace¬ 
ment  of  the  cyclonic  stominess  of  middle  latitudes  into  the 
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higher  latitudes.  The  poleward  displacement  of  the  descend¬ 
ing  branches  of  the  Ferrel  and  Hadley  cells  shifts  the  cur¬ 
rent  subtropical  high-pressure  belts  to  the  middle  latitudes; 
this  condition  yields  a  generally  settled,  mild,  storm-free 
climate  at  these  latitudes.  The  disappearance  of  the  easter¬ 
lies  at  high  latitudes  and  the  poleward  displacement  of  the 
westerlies  produce  a  relaxation  of  the  overall  circulation 
and  a  decrease  in  the  intensity  of  the  condensation  cycle. 

2.  For  glacial  conditions,  the  equatorward  displacement  of  the 
ascending  branches  of  the  polar  and  Ferrel  cells,  coupled 
with  a  marked  acceleration  of  the  condensation  cycle,  sus¬ 
tains  the  precipitation  of  higher  middle  latitudes  that 
builds  the  ice  sheets;  the  current  storm  tracks  of  middle 
latitude  roe  displaced  to  lower  latitudes,  and  they  exhibit 
increased  storminess  from  the  increased  intensity  of  the  cir¬ 
culation  in  the  Ferrel  cell.  The  acceleration  of  the  conden¬ 
sation  cycle  becomes  evident  from  the  precipitation  of  higher 
middle  latitudes,  the  extremely  pluvial  conditions  that  pre¬ 
vail  in  the  lower  middle  latitudes,  and  the  greatly  increased 
rainfalls  at  the  ITCZ. 

Finally,  the  climate,  as  characterized  by  the  mean  annual  surface 
temperature  T*(6),  for  the  current  climate  and  corresponding  tempera¬ 
tures  for  the  interglacial  and  glacial  conditions,  is  as  given  in 
Table  1  (Schell,  1961;  Fairbridge,  1961). 


TABLE  1.  T*(  0)  FOR  DIFFERENT  CLIMATES 


Latitude 

Present 

Interglacial 

Glacial 

80°-90°N 

-19°C 

0°C 

-70°C 

40°-60°N 

10 

13 

2 

0°-20°N 

26 

27 

23 

80°-90°S 

-49 

-30 

-110 
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Table  1  indicates  that  the  current  climate  is  closer  to  the  in¬ 
terglacial  climate  of  the  Pleistocene  Epoch.  The  table  also  shows  a 
rather  small  difference  between  the  interglacial  and  glacial  values 
of  T*  at  low  lat  tudes;  this  difference  increases  rapidly  at  the 
higher  latitudes  (Bell,  1953)  of  the  expanded  polar  cell. 

The  considerations  of  the  glacial-interglacial  conditions  indi¬ 
cate  that  the  cell  structure  of  the  current  general  circulation,  as 
induced  by  the  large-scale  transports  of  angular  momentum  and  sensible 
heat,  will  also  prevail  under  the  less  drastic  conditions  of  either 
advertent  or  inadvertent  climate  modification.  Also,  the  deduced 
general  circulation  for  the  extreme  glacial  and  interglacial  climates 
of  the  Pleistocene  Epoch  can  be  used  to  test  the  transitive  or  intran¬ 
sitive  character  of  a  nonlinear  model  that  relates  the  cell  structure 
of  the  general  circulation  to  the  net  heating  of  the  troposphere. 
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VI.  SAMPLE  OF  A  POTENTIAL  MECHANISM  FOR 
ADVERTENT  CLIMATE  MODIFICATION 

The  T*( 0)  values  at  the  middle  and  high  latitudes  for  the  inter¬ 
glacial  and  glacial  climates  of  the  Pleistocene  Epoch  indicate  that 
the  climate  of  the  earth's  surface  at  these  latitudes  is  very  sensi¬ 
tive  to  the  magnitude  of  the  surface  albedo,  A*.  Thus,  recalling  the 
third  and  fourth  basic  subquestions  identified  earlier,  it  is  of  in¬ 
terest  to  examine  briefly  if  man  has  achieved  the  technological  capa¬ 
bility  to  decrease  the  current  high  values  of  A*  at  the  high  northern 
latitudes. 

Figure  9  compares  the  magnitude  of  both  the  incoming  solar  radia¬ 
tion  and  its  absorption  by  the  earth's  surface  and  troposphere  at  low 
and  high  latitudes  (Budyko,  1966).  Even  though  the  solar  inclination 
is  low  at  high  latitudes,  the  length  of  solar  illumination  during  sum¬ 
mers  at  these  latitudes  is  such  that  the  magnitude  of  the  incoming 
solar  radiation  can  be  greater  than  that  at  low  latitudes  for  a  sig¬ 
nificant  fraction  of  each  year.  However,  because  of  the  relatively 
high  value  of  the  surface  albedo  at  the  high  latitudes,  the  absorbed 
short-wave  radiation  at  the  earth's  surface  at  these  latitudes  becomes 
significantly  lower  than  that  at  the  low  latitudes.  Thus,  a  decrease 
in  the  surface  albedo  accomplished  by  a  large-scale  removal  of  the 
snow  and  ice  at  the  high  latitudes  could  have  a  significant  effect  on 
the  structure  of  the  polar  and  Ferrel  cells  at  ta**  middle  and  high 
latitudes. 
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FIGURE  9.  Annual  Distribution  of  Incident  and  Absorbed  Solar  Radiation  at 
Equator  and  80°N  Latitude  (Budyko,  1966) 

The  feasibility  of  large-scale  removal  of  snow  and  ice  at  high 

latitudes  has  been  examined  theoretically  by  using  a  thermodynamic 

model  of  sea  ice  (Maykut,  1969),  since  sea  ice  is  a  dominant  feature 

7  2 

of  the  Arctic  Ocean,  which  covers  an  area  of  about  10  km  .  This  area 
is  about  80  percent  ice  during  the  winter  and  60  percent  in  the  late 
summer.  Because  of  the  large  horizontal  dimensions  of  the  surface  of 
the  Arctic  Ocean  relative  to  the  thickness  of  the  ice  and  snow  layers 
over  this  surface,  the  model  may  consider  an  infinite,  thin  slab  of 
ice  as  existing  over  the  water  and  an  infinite,  thin  slab  of  snow  as 
existing  over  the  ice.  Thus,  the  model  is  reduced  to  the  use  of  the 
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one-dimensional  heat  conduction  equation  for  both  the  snow  and  ice 
layers,  with  appropriate  boundary  conditions  for  each  layer.  The 
differential  equation  for  the  snow  layer  is  then  given  by 

(Pc)s  H  =  ks 

dz 

where  (k/pc)  is  the  thermal  diffusivity  for  the  snow,  T  the  tempera- 
ture  in  the  snow  layer,  z  the  vertical  coordinate,  and  t  the  time. 

The  differential  equation  for  the  sea  ice  must  account  for  the  varia¬ 
bility  of  the  thermal  diffusivity  k^/(  pc)^  from  the  inclusion  of  small 
pockets  of  salt  water  entrapped  when  sea  water  freezes.  Thus,  the 
product  of  the  density  and  specific  heat  (pc)^,  and  conductivity  k.  , 
for  sea  ice  may  be  expressed  as 


(  pc)i 


(po)i,o 


4100  Si(z) 
(T  -  273)2 


and 


0.28  Si(z) 
ki  "  ki,o  +  T "-—IT' 

where  the  subscript  i,o  denotes  pure  ice  and  a  known  ice  salinity 
profile.  Another  factor  that  must  be  considered  in  the  heat  conduc¬ 
tion  equation  for  the  snow  and  sea  ice  layers  is  the  absorption  of 
solar  radiation  per  unit  time,  per  unit  volume,  at  depth  z  in  the 
layers.  This  factor  may  be  expressed  by  using  the  extinction  coeffi¬ 
cient  h  as  hIq  exp  (-  hz),  where  IQ  denotes  the  amount  of  short-wave 
radiation  passing  through  the  surface  of  the  layers.  However,  for 
snow,  the  value  of  h  is  large  and  all  the  incoming  radiation  is  es- 
sentially  absorbed  in  the  first  few  centimeters  of  the  snow  layer. 

For  snow-free  conditions,  the  heat  conduction  equation  for  the  ice 
layer  is  then  given  by 


37 


(4-b) 


(  pc)t 


9T 

■St 


a2T 

j  +  h.Iq  exp  (-  k^z) 
Sz 


The  boundary  conditions  for  Eqs.  4-a  and  4-b  must  allow  for  the 
melting  of  snow  and  ice  at  their  interface  with  the  atmospheric  air, 
and  either  melting  or  accretion  of  ice  at  the  ice-water  interface. 

Thus,  the  boundary  conditions  are  as  follows: 

1.  At  the  snow-air  or  ice-air  interface,  the  net  energy  flux  is 
either  zero  when  T  <  273°K,  or  equal  to  the  heat  required  for 
melting  when  T  =  273°K. 

2.  At  the  snow-ice  interface  of  the  two  layers,  the  heat  flux  by 
conduction  in  the  snow  side  of  the  interface  must  be  equal  to 
that  in  the  ice  side  of  the  interface. 

3.  At  the  ice-water  interface,  the  net  heat  flux  from  conduction 
in  the  ice  side  of  the  interface  and  turbulence  from  the  ocean 
in  the  water  side  of  the  interface  must  equal  that  for  either 
melting  or  accretion  of  the  ice  at  this  interface. 

Since  the  parameters  at  the  interface  with  tie  atmospheric  air 
are  prescribed  as  a  function  of  time  (i.e.,  by  their  monthly  average 
values),  an  equilibrium  condition  is  eventually  reached  where  the  ice 
accretion  at  the  ice-water  interface  balances  the  ice  lost  through 
melting  at  the  ice-air  interface.  For  such  an  equilibrium  condition, 
the  net  mass  balance  for  the  year  is  zero,  and  both  the  mass  changes 
and  temperature  profiles  will  repeat  themselves  in  an  annual  cycle. 

The  numerical  integrations  of  Eqs.  4-a  and  4-b  yield  two  characteris¬ 
tics  of  the  ice  layer  that  are  of  interest  for  advertent  climate  modi¬ 
fication:  (1)  its  equilibrium  thickness  for  a  specified  monthly  at¬ 
mospheric  environment  during  the  year;  and  (2)  the  time  required  for 
a  change  from  one  equilibrium  condition  to  another  as  induced,  for 
example,  by  a  change  in  the  surface  albedo. 

Figure  10  shows  the  snow  and  ice  layers  over  the  water  for  the 
equilibrium  condition  deduced  from  Eqs.  4-a  and  4-b  for  the  current 
atmospheric  and  oceanic  environments.  The  specified  snow  layer  is 
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based  on  experimental  data;  after  the  freeze-up  in  late  August,  it  is 
assumed  to  grow  in  thickness  to  30  cm  by  the  end  of  October  and  to  40 
cm  during  May.  The  figure  shows  no  snow  during  the  melt  season.  The 
equilibrium  thickness  of  the  ice  is  determined  by  using  arbitrary  ini¬ 
tial  conditions  for  the  temperature  and  ice  thickness  in  the  integra¬ 
tion  of  the  equations.  This  figure  shows  that  the  model  predicts  an 
average  equilibrium  ice  thickness  for  the  current  environments  of 
about  3  m,  a  value  that  agrees  rather  well  with  measurements  from 
drifting  stations  on  pack  ice.  There  is  also  good  agreement  between 
the  monthly  calculated  and  observed  surface  temperatures,  which  yield 
annual  values  of  -18.2°C  and  -18.6°C,  respectively.* 


FIGURE  10.  Annual  Ice  Equilibrium  Thickness  and  Temperature  Reid  as 
Determined  from  One-Dimensional  Transient  Thermodynamic 
Model  of  Sea  Ice  (Maykut,  1969) 


Because  of  the  relatively  small  thickness  of  the  ice  and  snow  layers, 
the  argument  that  the  melting  of  the  ice  and  snow  at  the  poles  would 
flood  cities  does  not  apply  to  the  Arctic  Ocean. 
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The  model  can  be  used  to  investigate  systematically  the  relative 

effect  of  changes  in  each  variable  defining  the  atmospheric  and  oceanic 

environments  on  the  equilibrium  ice  thickness.  Figure  11  shows,  for 

example,  the  effect  bn  the  equilibrium  ice  thickness  of  changes  in 

only  the  oceanic  environment,  vdiich  is  defined  by  a  constant  heat  flux 

from  the  ocean  throughout  the  year  at  the  ice -water  interface  and 

equals  its  annual  average.  It  also  shows  the  maximum,  average,  and 

minimum  equilibrium  ice  thickness  (as  exhibited,  for  example,  in  Fig. 

10)  as  the  magnitude  of  the  oceanic  heat  flux  is  varied  parametrically 

between  0  and  4.5  kcal/cm  year.  As  indicated  by  the  figure,  the 

equilibrium  ice  thickness  would  almost  double  its  current  value  if  the 

oceanic  heat  flux  would  vanish;  the  numerical  results  also  indicate 

that  the  ice  would  vanish  when  the  oceanic  heat  flux  is  increased  to 

2 

about  5  to  6  kcal/cm  year.  The  relative  results  in  Fig.  11  are  use¬ 
ful  when  attempting  to  evaluate  the  merit  Of  inducing  warmer  water 
into  the  Arctic  for  the  purpose  of  large-scale  ice  removal.  However, 
systematic  investigations  with  this  model  indicate  that  the  equilibrium 
ice  thickness  is  most  sensitive  to  the  magnitude  of  the  surface  albedo. 
These  investigations  indicate  the  following  results:  (1)  if  the  cur¬ 
rent  atmospheric  and  oceanic  environments  are  kept  unchanged  except 
for  a  10  percent  reduction  in  the  surface  albedo  during  June  1  to  the 
onset  of  the  snow  only  (Fig.  10),  the  equilibrium  ice  thickness  is 
reduced  to  about  1  m  or  1/3  of  its  current  value;  and  (2)  if  the  sur¬ 
face  albedo  is  reduced  by  20  percent  instead  of  the  10  percent  in  the 
previous  case,  the  ice  vanishes  during  the  summer  of  the  third  year. 
Since  measurements  from  high  towers,  and  calculations  using  satellite 
data,  indicate  a  large-scale  average  of  about  0.5  for  the  summer  sur¬ 
face  albedo,  these  relative  results  suggest  that  the  use,  over  the 
Arctic  surface,  of  a  suitable  material  (i.e.,  lighter  than  water, 
slowly  soluble  in  water,  nontoxic,  etc.),  dark  enough  to  yield  a 
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large-scale  summer  albedo  of  about  0.4,  would  eliminate  the  ice  layer 
over  the  Arctic  within  very  few  years.* 


FIGURE  11 .  Equilibrium  Ice  Thickness  as  a  Function  of  the  Oceanic  Heat  Flux 
to  the  Sea  Ice  Layer  as  Determined  Parametrically  from  Thermo¬ 
dynamic  Model  of  Sea  Ice  (Maykut,  1969) 


Since  the  model  ignores  the  effect  of  melt  ponds  and  leads  (which 
have  a  lower  surface  albedo  than  that  of  ice),  the  experimental 
large-scale  averages  of  the  surface  albedo  are  somewhat  lower  than 
those  for  melting  ice  used  in  the  model. 
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VI I. 


DESCRIPTION  OF  THE  GENERAL  CIRCULATION  FOR  CLIMATE 


The  large-scale  removal  of  the  sea  ice  over  the  Arctic  Ocean 
would  have  an  impact  on  the  net  heating  of  the  troposphere  at  high 
latitudes,  through  changes  in  each  of  the  terms  in  Eq.  1.  The  change 
in  the  net  heating  of  the  troposphere  at  high  latitudes  would  modify 
the  current  polar  cell,  which  in  turn  would  affect  the  atmospheric 
environment  at  these  latitudes.  The  description  of  the  coupling  be¬ 
tween  the  long-term  general  circulation  and  net  heating  of  the  tropo¬ 
sphere  becomes,  then,  essential  in  establishing  both  the  magnitude  and 
the  transitive  or  intransitive  nature  of  the  changes  of  the  climate  at 
the  earth’s  surface,  from  advertent  modifications  of  parameters  such 
as  the  surface  albedo. 

The  statistics  of  the  general  circulation  for  climate  can  be  de¬ 
termined  from  the  statistics  of  the  net  heating  of  the  troposphere  by 
two  approaches: 

1.  Emphasize  the  rather  straightforward  numerical  integrations 
of  the  equations  of  motion  (e.g.,  Holloway  and  Manabe,  1971; 
Kasahara  and  Washington,  1971;  Arakawa,  1966;  Mintz,  1964). 
This  approach  is  an  adaptation  of  the  numerical  methods  de¬ 
veloped  for  long-range  weather  forecasting  to  climate,  and 
it  requires  the  use  of  very  large  electronic  computers  to 
perform  very  long  numerical  experiments.  It  therefore  yields 
a  very  detailed  description  of  atmospheric  phenomena,  which 
are  derived  from  difference  equations  approximating  the  dif¬ 
ferential  equations  of  motion.  The  limit  on  the  capability 
of  a  given  computer  introduces  an  upper  limit  to  the  wave 
number  that  may  be  considered  by  the  difference  equations. 
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2.  Emphasize  the  parameterization  of  physical  phenomena  in  the 
numerical  integration  of  the  equations  of  motion  (e.g., 
Smagorinsky,  1964;  Kurihara,  1970;  Saltzman  and  Vemokar, 
1971).  This  approach  requires  maximum  physical  insight  into 
the  key  processes  controlling  tropospheric  phenomena,  e.g., 
the  formulation  of  the  countergradient  transport  of  angular 
momentum  as  well  as  the  transports  of  sensible  heat  and  hu¬ 
midity  by  the  large-scale  two-dimensional  turbulence. 

As  compared  with  the  former  approach,  the  latter  approach  can 
carry  out  the  numerical  integrations  of  the  equations  of  motion 
through  the  use  of  smaller  electronic  computers,  and  the  essential 
long-term  statistics  of  the  general  circulation  can  be  derived 
without  the  need  of  the  very  long  time  integrations  that  are  required 
in  the  former  approach.  For  these  reasons  the  effort  in  this  paper 
is  concerned  exclusively  with  tills  latter  approach. 

The  understanding  of  the  three-cell  structure  of  the  general  cir¬ 
culation  for  the  current  climate  was  achieved  only  rather  recently, 
and  it  was  derived  from  the  experimental  evaluation  of  the  large-scale 
transport  of  angular  momentum  at  the  upper  levels  of  the  troposphere 
(Lorenz,  1967).  These  investigations  indicated  that  the  streamlines 
(♦*)  of  the  flow  in  the  troposphere  at  middle  latitudes  are  as  sketched 
in  Fig.  12,  i.e.,  the  motion  is  defined  by  a  westerly  motion  on  which 
troughs  in  the  general  vicinity  of  cyclonic  disturbances  are  super¬ 
posed  (Starr,  1948).  The  important  feature  of  these  streamlines  is 
their  departure  from  a  symmetrical  sinusoidal  form  and  the  associated 
NE-SW  tilt  from  the  meridian  (a)  of  the  trough  line.  This  departure 
from  sinusoidal  form  or  tilt  of  the  trough  line  is  a  function  of 
latitude,  i.e.,  the  tilt  is  small  at  higher  latitudes  near  the  north¬ 
ern  border  of  the  surface  westerly  belt  (Fig.  5)  vhere  the  transport 
of  angular  momentum  is  small,  and  it  becomes  more  pronounced  at  lower 
latitudes  vhere  the  transport  of  angular  momentum  is  large.  Figure  12 
shows  also  the  temperature  disturbances  ($*),  vhich  can  be  defined  by 
introducing  a  longitudinal  phase  lag  (6)  between  the  temperature  and 
flow  disturbances  (Smagorinsky,  1964).  The  analytical  representations 
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of  these  disturbances  are  included  in  Fig.  12,  and  they  are  defined  by 
the  tilt,  the  thermal  phase  lag,  a  zonal  wave  number  appropriate  for 
climate,  and  the  disturbance's  respective  amplitudes  A  and  B;  thus, 
these  formulations  possess  the  minimum  degrees  of  freedom  necessary  to 
effect  the  poleward  transfer  of  both  angular  momentum  and  sensible 
heat.  A  brief  but  coherent  description  of  the  basic  steps  in  the  cal¬ 
culation  of  these  statistics  (e.g.,  the  annual,  zonal  averages  of  the 
angular  momentum  and  sensible  heat  transports)  for  climate  from  the 
statistics  of  the  net  heating  rotes  of  the  troposphere  (Eq.  1)  is 
given  below. 


t*  -  A  tin  k(X  -»e) 

S  *“  8  iln  M  X  -  98  *6) 


FIGURE  12.  Schematic  Diagram  of  a  Baroclinic  Disturbance 
(Starr,  1948;  Smogorinsky,  1964) 


The  basic  conservation  equations  for  the  baroclinic  flow  in  the 
troposphere  are  as  follows  (e.g.,  Lorenz,  1967): 
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(S-a) 


dV/dt  = 

dT/dt  = 

d  p_1/dt  = 

P  = 

These  equations  represent, 
or  momentum  equation,  the  first  law  of  thermodynamics  or  energy  equa¬ 
tion,  the  conservation  of  mass  or  continuity  equation,  and  the  equation 
of  state  for  air.  The  left-hand  side  of  Eq.  5-a  denotes  acceleration 
of  the  three-dimensional  wind  d;  the  terms  in  the  right-hand  side  are, 
respectively,  the  Coriolis  acceleration,  the  pressure  gradient,  the 
earth's  gravity,  and  the  friction  force  per  unit  mass.  In  Eq.  5-b,  the 
symbol  cp/cy  denotes  the  ratio  of  the  specific  heat  of  air  at  constant 
pressure  to  the  specific  heat  at  constant  volume,  and  Q  denotes  net 
heating  per  unit  mass.  Using  the  conditions  of  hydrostatic  equilibrium 
or  no  vertical  wind,  Eqs.  5-a  through  5-d  can  be  put  in  the  so-called 
primitive  form,  where  the  wind  field  is  denoted  by  a  two-dimensional 
vector  (d),  and  the  independent  variable  z  in  the  radial  direction 
from  the  earth's  center  is  replaced  by  the  pressure  p.  The  primitive 
equations  are  as  follows  (e.g.,  Lorenz,  1967): 


•20  xV-~7p  +  g  +  F 


4H 


T  7  •  V  +  Q/c 


(5-b) 


I 

-  V  •  V 
p 


( 5-c) 


pRT 


(5-d) 


respectively,  Newton's  second  law  of  motion 


dU/dt  =  -f  k  x  U  ■  g  Vz  +  F 

(5-e) 

dT/dt  =  RT<i>/CpP  Q/Cp 

(  5— f  ) 

^  •  U  +  du)/dp  =  0 

(5-g) 

dz/dp  =  -RT/gp 

( 5-h) 
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where  U  =  iu  +  jv,  with  i  and  j  denoting  unit  vectors  in  the  eastward 
and  northward  directions,  f  is  the  Coriolis  parameter  20  sin  6,  £  the 
unit  vector  i.  the  radial  direction,  9  is  now  a  horizontal  differen- 
t'al  ope^tor  holding  the  pressure  constant,  and  uu  =  dp/dt.  Also,  to 
make  the  kinetic  energy  of  the  horizontal  motion  consistent  with  the 
hydrostatic  approximation,  dr/dt  =  0  or  r  must  be  replaced  by  the 
earth’s  radius  a  (Smagorinsky ,  1958).  Equation  5-e  describes  the 
horizontal  motion  and  it  yields  two  scalar  components,  one  for  the 
zonal  flow  and  the  other  for  the  meridional  flow;  Eq.  5-f  is  the  en¬ 
ergy  equation,  Eq.  5-g  the  continuity  equation,  and  Eq.  5-h  is  the 
hydrostatic  condition  using  Eq.  5-d.  The  continuity  equation  is  put 
in  the  form  5-g  from  the  ’ise  of  both  the  hydrostatic  equation  to  elim¬ 
inate  the  density  p,  and  the  transformation  from  vertical  (z)  to  the 
pressure  (p)  coordinate  (e.g.,  Thompson,  1961).  In  the  spherical  co¬ 
ordinates  of  the  earth  (X,  0),  the  two  scalar  equations  for  the  zonal 
and  meridional  motion  are,  respectively,  as  follows  (e.g.,  Lorenz, 
1967): 


(Si) 

a£  *  •  T  tan  8  ■  fu  ■  ?  H  *  Fe  (5*)) 

where  u  =  a  cos  0  dX/dt,  v  =  a  d0/dt,  and  0  denotes  the  geopotential 
gz.  The  total  or  individual  time  derivatives  in  Eqs.  5-i,  5-j,  and 
5-f,  say  dX/dt,  are  related  to  their  local  time  derivatives  (dX/dt) 
by  the  equation 


dX  _  ax 
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The  primitive  Eqs.  5-e  through  5-h  can  be  applied  to  a  two-level 
model  of  the  atmosphere.  This  model  represents  the  flow  pattern  in 
the  upper  half  of  the  atmosphere  by  the  flow  pattern  at  the  250-mb 
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level,  and  the  flow  pattern  in  the  lower  half  of  the  atmosphere  is 
represented  by  the  flow  pattern  at  the  750-mb  level  (Chamey  and 
Phillips,  1953;  Phillips  1956),  Figure  13  shows  these  two  levels, 
together  with  the  dependent  variables  that  are  applicable  at  these 
levels.  The  momentum  and  continuity  equations  are  applied  at  the 
2  °“mb  and  750-mb  levels,  while  the  energy  equation  is  applied  at  the 
500-mb  level  (as  implied  previously  in  Fig.  1)  together  with  the  hy¬ 
drostatic  equation  5-h.  Figure  13  also  shows  the  boundary  conditions 
at  the  top  and  bottom  of  the  atmosphere,  i.e.,  «^  =  u>4  =  0,  which  are 
used  to  filter  gravity  waves. 
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FIGURE  13.  Two- Level  Model  of  Borocllnic  Flowi  (Phillip*,  1956) 


The  large-scale  transports  of  angular  momentum  and  sensible 
heat  were  calculated  by  using:  (1)  the  primitive  equations,  5-i,  5-f, 
5-g,  and  5-h  but  not  5-j;  (2)  the  two-level  model  of  the  atmosphere 
(Fig.  13);  and  (3)  the  mechanisms  of  the  tilt  and  phase  lag  described 
in  Fig.  12.  The  calculation  of  these  transports  was  reduced  to  the 
solution  of  a  two-parameter  eigenvalue  problem,  which  involved  the 
numerical  integration  of  a  nonlinear,  ordinary,  second-order  differ¬ 
ential  equation  (Smagorinsky,  1958,  1963,  1964).  Although  a  unified 
and  somevfcat  detailed  description  of  this  procedure  is  included  in 
Appendix  A  to  clarify  all  the  approximations  used  in  this  model,  a 
broad  outline  of  steps  (1)  through  (3)  is  given  by  the  following  two 
results: 
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A.  TRANSPORTS  AS  A  FUNCTION  OF  POLEWARD  HEAT  FLUX 


The  use  of  the  zonal  momentum  Equation  5-i,  together  with  the 
continuity  Equation  5-g  at  the  250-mb  and  750-mb  levels,  yields  two 
zonal  equations,  which  after  taking  a  zonal  average  of  each  defined 
by 


[X]  * 
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£na 


X  dx 


can  be  written  as  (Appendix  A): 
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Equations  5-k  and  5-t  are  written  in  the  Mercator  (x,y)  coordinates  in¬ 
stead  of  the  spherical  (X,0)  coordinates  by  using  dx  =  a  dX  and  dy  = 
a  dB/cos  0.  The  zonal  and  meridional  winds  in  the  Mercator  plane  are 
related  to  their  counterparts  in  the  spherical  coordinates  through  the 
relations  um  -  urn  and  vm  =  vm,  where  m  =  1/cos  0.  The  notation  X  and 
X  denote,  respectively,  the  sum  and  the  difference  of  the  given  param¬ 
eters  at  levels  one  and  three  (Fig.  13),  i.e.,  X  =  X^  +  Xj  and  X  = 

Xx  -  X^;  Mm  denotes  a  quantity  proportional  to  the  eddy  flux  of  an¬ 
gular  momentum,  i.e., 


Finally,  the  symbols  C  and  D  in  Eq.  5-t  denote,  respectively,  the  rela¬ 
tive  vorticity  and  the  horizontal  divergency,  which  are  given  by 
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C  =  dvm/dx  “  ro2  d(um/m2)/dy 


Similarly,  the  use  of  the  energy  Equation  5-f  and  the  hydrostatic 
Equation  5-h  at  the  500-mb  level  (i.e.,  Fig.  1)  yields  the  relation 
(Appendix  A): 
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where  Hm  denotes  a  quantity  proportional  to  the  eddy  fl  x  of  sensible 
heat,  i.e., 


H»  = 


~ 7“ 


2 

and  y 


is  defined  by 


Y2  =  (l/2)p2(d*/dp)(d  In  Qp  /dp) 


where  0p  is  the  potential  temperature,  i.e., 


R/c 

9p  =  (P4/P) 


P  . 
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Introducing  the  condition  of  quasi-equilibrium  in  Eqs.  5-k,  5-4,  and 
5-m,  the  local  time  operator  vanishes,  i.e.,  d/dt  =  0.  Neglecting 
the  first  term  in  the  right-hand  side  of  Eq.  5-k,  the  Q  and  [D*u^]/m2 
in  Eq.  5-4,  and  assuming  that  the  frictional  forces  in  Eqs.  5-k  and 
5-4  are  due  only  to  the  stresses  at  the  earth’s  surface  ( t^)  or 
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[F]  -  -[r ]  -  g[T4]/p,  this  system  of  equations  reduces  to  the  follow¬ 
ing  form: 


"'4  IjT  =  [tx4^ 
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Since  y  is  a  constant  related  to  the  static  stability,  integration 
of  Eq.  5-p  yields 


where 


=  qf(y) 


qf(y)  = 


dy 


Using  the  approximation  M  Mm  ±  in  Eq.  5-n,  and  using  Eq.  5-o 

to  eliminate  [v  ]: 
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2m2  3M, 
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Thus,  the  equations  of  motion  yield  Eq.  5-q  for  a  relationship  involv¬ 
ing  the  transport  of  sensible  heat,  Hm,  the  latitudinal  gradient  of 
the  transport  of  angular  momentum,  Mm,  the  annual  mean  energy  flux, 
qf ,  as  obtained  from  the  net  heating  rate  of  the  atmospheric  column 
(i.e.,  [Q2J  as  given  by  Eq.  1),  and  one  eigenvalue  parameter  y2. 
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B.  TRANSPORTS  AS  A  FUNCTION  OF  WAVE  PARAMETERS 

The  use  of  the  mechanism  of  the  tilt  and  phase  lag  described 
in  Fig.  12,  together  with  the  condition  of  geostrosphic  equilibrium, 
leads  to  another  relation  between  the  transports  of  sensible  heat 
and  angular  momentum.  The  geostrophic  equilibrium  is  a  feature  of 
the  circulation  in  middle  and  higher  latitudes  that  is  almost  as 
prominent  as  hydrostatic  equilibrium;  it  is  given  by  the  approximate 
balance  between  the  Coriolis  force  and  the  horizontal  pressure  grad¬ 
ient  in  Eq.  5-e,  a  condition  that  yields  -fic  x  U  =  ^0,  or  equiva¬ 
lently,  U  =  (1/f )k  x  ^0.  The  geostrophic  equilibrium  allows  the 
zonal  and  meridional  eddy  winds  to  be  expressed  in  terns  of  the  tem¬ 
perature  disturbance  in  Fig.  12  as  follows: 


,?*  _  90*  _  m  90* 

m  T  "5y - aF  W 


(5-r) 


90*  _  m2  90k 
m  T~  9x  aF  bx 


(5-s) 


The  continuity  Equation  5-g,  with  the  boundary  conditions  on  (Ug  and 
(Fig.  13),  yields  a  relationship  between  the  zonal  and  meridional 
eddy  winds  and  the  flow  disturbances  in  Fig.  12  as  follows: 


m  9\|r* 
a  Te 


9**  _  m2  9** 
Sx  a  'd\ 


whe re,  as  indicated  in  Fig.  12, 


t*  =  A  sin  71 

0*  =  B  sin  (7)  +  k6) 


(5-t) 


(5-u) 


(5-v) 
(  5-w) 
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with  T|  =  k(  X  -  aQ).  The  second  relationship  between  the  transports  of 
sensible  heat  and  angular  momentum  is  then  obtained  from  differentia¬ 
tion  of  the  expression  for  H  ,  use  of  the  definition  of  M  ,  and  the 

m  m 

use  of  the  relations  5-r,  5-u,  and  5-w;  the  results  can  be  written  as 
follows  (Appendix  A): 


M 


m 


K 


dHm 


where  S  is  a  second  eigenvalue  parameter  given  by 


(5-x) 


a  tan  k6 

— — 


This  definition  of  S  stems  from  the  evaluation  of  the  terms  [u*^]/m4 
and  Hm  from  the  relations  5-r,  5-u,  and  5-w,  i.e., 


„  _  kAB 

Hm  "  3a“ 


[sin  ( T)  +  k 6)  cos 


’ll  =  Sir 


sin  k6 


(5-y) 


and 


CW 


[cos  („  +  kW  cos  „]  =  k  k» 

2a  rm 


m 


"2. 
a  fm 


This  latter  expression  is  derived  by  assuming  that  B,  or,  and  6  in 
Eq.  5-w  are  nearly  constant;  thus,  Eq.  5-y  indicates  then  that  Hm  is 
proportional  to  A.  The  combination  of  these  two  latter  expressions 
yields  the  definition  of  S,  since  [u*vj]/m4  can  be  written  as  follows: 


[u*9*  ]  H 
u  m  m  J  m 

3  fmS 

m 
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The  combination  of  the  two  foregoing  relations  5-q  and  5-x  for 
the  transports  of  sensible  heat  and  angular  momentum  yields: 


If  S  is  assumed  to  be  constant  given  by  a  proper  average,  Fq.  5-z  is 

a  second-order,  ordinary  differential  equation  with  two  eigenvalue 

2 

parameters  given  by  the  constants  y  and  S;  these  parameters  are  fixed 
from  the  use  of  four  boundary  conditions,  i.e.,  Hm(0,y)  =  0  and 
dHm/dy  =  0  at  y  =  0,  y  =  Y,  where  Y  denotes  a  latitude  near  the  poles. 

Either  Eqs.  5-q  and  5-x  or  5-z  and  5-x  yield  the  large-scale 
transports  of  sensible  heat  and  angular  momentum  as  a  function  of  the 
poleward  heat  flux  (q^)  induced  by  the  net  heating  rates  of  the  at¬ 
mospheric  column  (i.e.,  Eq.  1).  The  identification  of  q^  as  a  pole- 
ward  transport  of  energy  is  formally  derived  by  integrating  the  prod¬ 
uct  of  a  mass  element  and  the  individual  time  derivative  of  any  de¬ 
pendent  variable  X  over  the  volume  of  the  region  north  of  a  latitude 
Gp  and  then  taking  the  time  and  zonal  averages  in  such  an  integra¬ 
tion.  Thus,  the  use  of  the  definition  for  the  individual  time  deriva¬ 
tive  of  X,  the  continuity  Equation  5-g  in  spherical  coordinates  (\, 

-1  2 

0),  the  mass  element  g  a  cos  0  d\  d0  dp,  and  taking  the  time  and 
zonal  averages  yield  the  following  expression  (e.g.,  Lorenz,  1967): 


2na  cos  0-^CXvjg"1  dp  = 


cos0[dX/dt]g_±  d0  dp 


2  2 

Letting  X  denote  the  energy  per  unit  mass  e(cal/gm  or  m  /sec  ),  Q  the 
net  heating  rate  de/dt,  and  Q2  =J  Q  d(p/p^),  Eq.  6  yields 
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/ 


—  "I 


2m  cos  e^ev]  g  dp 


P^2 

g 


cos  8  de 


Using  the  foregoing  definition  of  q^,  the  integral  in  the  right-hand 
side  of  this  expression  may  be  written  as 


f 


-1 


2rra  cos  [ev]  g  dp  - 


- 2 ;ft/5~~7ag  |Mn/2)  - 


or,  since  q^(TT/2)  -  0> 


qf(61)  =  I  [ev]  cos  0X  d(p/p4) 

Jg  P 

which  relates  qf( 0)  to  the  average  poleward  transport  of  energy  [ev] 
from  the  atmospheric  column  at  a  latitude  0^.  Figure  14  shows  the 
function  qf  corresponding  to  the  processes  indicated  in  Fig.  2  for 
the  current  climate.  Figure  14  also  shows  schematically  the  modifi¬ 
cation  of  qf  at  the  higher  latitudes  by  the  advertent  removal  of  the 
ice  from  the  Arctic  Ocean;  this  modification  would,  of  course,  result 
from  the  corresponding  changes  in  each  of  the  terms  in  Eq.  1  from  the 
advertent  change  of  the  surface  albedo  A*  at  these  higher  latitudes. 
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CUMFNT  CLIMATE 


NORTH  LATITUDE,  d«g 


FIGURE  14.  Poleward  Heat  Flux  Required  by  the  Net  Heating  Rate,  d aQJq 
(Smag  or  insky,  1963)  ^  * 
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VIII .  EVALUATION  OF  THE  THEORY  OF  THE  GENERAL  CIRCULATION 
FOR  CLIMATE  BASED  ON  A  CONSTANT  S  WAVE  PARAMETER 


The  use  of  Eqs.  5-z  and  5-x  with  a  constant  S  wave  parameter 
yields  predictions  for  the  two  important  large-scale  transports  of 
sensible  heat  and  angular  momentum  as  a  function  of  latitude,  i.e., 
the  transports  that  control  the  meridional  cell  structure  of  the  gen¬ 
eral  circulation  (Fig.  5).  The  results  from  this  model  can  be  evalu¬ 
ated  from: 

1.  A  comparison  of  the  predicted  transports  with  corresponding 
experimental  values  obtained  from  the  reduction  of  data  con¬ 
cerning  the  statistics  of  the  general  circulation. 

2.  Theoretical  considerations  concerning  the  sensitivity  of  the 
predicted  values  of  these  transports  to  the  characteristics 
of  the  S  wave  parameter,  vdiich  is  determined  from  boundary 
conditions  at  the  north  pole. 

3.  A  comparison  of  the  semi -experimental  behavior  as  a  function 
of  latitude  of  the  amplitude  and  wave  parameters  of  the 
baroclinic  disturbances  with  their  assumed  behavior  in  the 
model. 

Equations  5-q  and  5-x  yield  power  series  solutions  for  near  the 
equator,  which  are 
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where 


«f 


e-*o 


e 


and  power  series  solutions  for  near  the  poles: 


vfriere 


C  and  r  are  constants,  and  r  Js  1  to  satisfy  the  boundary  condition 
vm  The  Power  series  solutions  for  near  the  poles  yield 

values  for  both  H^Y)  and  dHm/dy  as  y  -*  Y,  v^iich  are  required  for  the 
determination  of  the  eigenvalue  parameters  and  S  from  the  numeri¬ 
cal  solutions  of  Eqs.  5-q  and  5-x  in  the  range  0  s  y  s  Y,  or  equiva¬ 
lently  0  s  e  £  tt/2. 

Figure  15  shows  the  large-scale  transport  of  sensible  heat,  Hm, 
derived  from  the  numerical  integrations  of  Eqs.  5-q  and  5-x  and  using 
the  function  q^  in  Fig.  14  for  the  current  climate.  Table  2  shows 
experimental  annual  values  for  the  northward  transports  of  sensible 
heat  (CpT)  by  the  transient  [v'T']  and  stationary  [v*T*]  eddies  as 
obtained  from  the  most  recent  reduction  of  available  data  (Oort  and 
Rasmusson,  1971).  Using  the  definition  of  Hm,  the  transformation 
vm  =  vm,  the  relation  =  RT£  and  9*  =  2v£,  give  Hm  =  { [v^T^]  + 

[v*T* cos  Since  Q2  is  defined  by  d(p/p4)  in  Eq.  1,  the 


58 


NORTH  LATITUDE,  d*s 


FIGURE  15.  Comparison  Between  Theory  (Smagori nsky,  1969)  and  Experimental 
Annual  Data  (Oort-Rasmusson,  1971)  for  the  Transient  Plus 
Standing  Eddy  Transport  of  Heat 


TABLE  2.  NORTHWARD  ANNUAL  TRANSPORT  OF  SENSIBLE  HEAT 
BY  TRANSIENT  EDDIES  [v  'T  ']  AND  STATIONARY 
EDDIES  [v*T* ],  m°C/sec,  AT  LATITUDES 
0  s  e  s  75°N 

Transient  Eddies 


Stationary  Eddies 


transports  within  the  bracket  ar«  also  given  by  the  mean  vertical 
values  of  the  atmospheric  column,  i.e.,  [v^T^ ]  =  JgCv'T']  d(p/p4)  as 
well  as  [v£T|]  =  d(p/p4).  The  experimental  values  of  Hm  ob¬ 

tained  in  this  manner  from  the  mean  values  in  Table  2  are  also  shown 
in  Fig.  15;  this  figure  indicates  a  good  agreement  between  the  theo¬ 
retical  and  experimental  values  of  H  .  Figure  16  shows  similar  data 
for  the  large-scale  transport  of  angular  momentum  at-  the  250-mb  levei., 

i.e.,  ..  The  experimental  values  of  M  ,  are  obtained  by 

m  m,l  r  m,l  7 

using  the  definition  of  M  .  the  transformations  u  =  urn  and  v  =  vm: 

—r—  _m_  2  m  m  ’ 

thus,  M  .  =  { [u V ']  +  [u*v* jJcos^G,  where  the  transports  are  shown  in 

Table  3.  Figure  16  indicates  that  the  theory  underestimates  the 
transport  of  angular  momentum  at  midlatitudes  by  as  much  as  about  40 
percent.  Figure  16  shows  a  good  agreement  for  the  latitudes  of  the 
largest  positive  and  negative  absolute  magnitudes  of  the  angular  mo¬ 
mentum  transports,  or  the  latitudinal  widths  of  the  Hadley  and  Ferrel 

cells,  i.e.,  where  dM  ,/d0  =  0.  This  latter  condition  is  obtained 

m  ,  l 

from  Eqs.  5-n  and  5-o,  which  yield 

[vm]  =  (2m3/af)  •  dMm  /d0 

since  the  continuity  equation,  together  with  the  boundary  conditions 

for  the  meridional  wind,  vm  =  0  at  the  equator  and  the  poles  at  both 

levels  (Fig.  13),  gives  [9  ]  =  0,  i.e.,  [v  ]  =  -[v  ].  Then 

m  *1*1  ***  3 

[v  ]  =  (m3/af)  •  dM  /d0 

m  j 

Thus,  the  result  [v  n]  =  -[v  ,]  indicates  that  the  meridional  cir- 

culations  in  the  upper  and  lower  branches  of  each  cell  (Fig.  5)  are 
equal  but  opposite  in  direction,  vshile  the  two  conditions  dMm  ^/d0  =  0 
in  Fig.  16  yield  the  three -cell  structure  of  the  current  long-term 
general  circulation. 
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TRANSPORT  OF  ANGULAR  MOMENTUM, 


[ 
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TABLE  3.  NORTHWARD  ANNUAL  TRANSPORT  OF  WESTERLY  MOMENTUM 
BY  TRANSIENT  EDDIES  [v  'u  ']  AND  STATIONARY 
EDDIES  [v*u*],  m2/sec2,  AT  LATITUDES 
0  s  0  s  75°N 

Transient  Eddies 
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0 

5 

10 

15 

20 

25 

30 

35 

40 
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SO 

m 

60 

65 

70 

75 

SO 
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D 
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ta 
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38.0 
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-0.1 

-3.9 

-5.3 
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300 
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m 
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m 
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Stationary  Eddies 
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If  Eqs.  5-x  and  5-z  are  multiplied  by  an  arbitrary  constant,  both 
the  transports  of  angular  momentum  and  sensible  heat  are  directly  pro- 

a 

portional  to  q^.  Thus,  the  magnitude  of  the  parameters  y  and  S  re- 
mains  fixed  for  a  given  shape  of  the  q^  curve,  and  the  discrepancy  be¬ 
tween  theory  and  experiment  for  the  angular  momentum  transport  in  Fig. 
16  could  not  be  corrected  by  an  increase  in  the  magnitude  of  qf  with¬ 
out  degrading  the  agreement  for  the  transport  of  sensible  heat  in  the 
previous  figure. 

The  foregoing  constant  S  theory  of  the  general  circulation  for 
climate  is  successful  in  the  accurate  prediction  of  the  large-scale 
transport  of  sensible  heat,  and  in  the  determination  with  a  lower  de¬ 
gree  of  accuracy  of  the  countergradient  transport  of  angular  momentum 
for  the  upper  levels  of  the  troposphere  at  middle  and  higher  latitudes 
where  d[u]/d0  <  0.  However,  a  more  decisive  test  of  this  theory  can 
be  made  from  both  theoretical  considerations  and  the  use  of  experi¬ 
mental  data  for  other  statistics  of  the  long-term  general  circulation. 

To  test  the  sensitivity  of  the  transports'  magnitude  to  the 

boundary  conditions  that  are  used  to  fix  the  constant  S,  an  arbitrary 

variable  S  (0),  defined  as  S  sin  20,  can  be  used  in  either  Eq.  5-x 
p  max 

or  5-z.  The  expression  for  S  ( 0)  has  a  maximum  value  S  at  45-deg 

p  max 

latitude  and  it  becomes  zero  at  both  the  equator  and  the  north  pole. 

The  boundary  conditions  for  H  and  dH  /d0  as  0-60  deg  are  now  used  to 

2  ni  m 

determine  the  constants  y  and  S  .  Numerical  integrations  of  Eqs. 

max 

5-q  and  5-x,  using  Sp( 0)  instead  of  a  constant  S,  indicate  that  the 
transports  are  insensitive  to  Sp(0);  i.e,,  the  theoretical  results 
shown  in  Figs.  15  and  16  apply  essentially  to  both  Sp( 0)  and  a  con¬ 
stant  S .  The  use  of  Sp( 0)  yields  maximum  values  for  Hm  and  ^  that 

are,  respectively,  only  about  2  percent  and  4  percent  higher  than  the 

2 

corresponding  values  for  a  constant  S;  the  value  of  y  for  Sp( 0)  is 

only  2  percent  higher  than  that  for  a  constant  S,  while  the  value  of 

S  „  is  about- 5  percent  lower  than  the  constant  S  value.  Thus,  this 
max 

particular  change  in  S  implies  that  the  transports  at  middle  latitudes 
can  be  insensitive  to  the  boundary  condition  dHm/de  at  the  pole,  and 
that  only  the  tilt  and  phase  lags  at  middle  latitudes  are  important. 
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3. 


are  *"  th*  «“  *  *•  — -  s  theory 

1.  It  cannot  verify  the  assumptions  used  In  the  derivation  of 
ssentrany  the  parameter  S  of  the  behavior  of  the  parame- 

*  and  6’  whlch  are  held  constant.  These  latter  as¬ 
sumptions  become  Important  for  the  verification  of  a  consis- 
tcncy  between  dH  /de  as  derived  from  Ea  5-v  ->nH  h, 
shown  m  Fig.  15.  5  y  Md  Che  results 

The  use  of  the  theory  for  very  low  latitudes,  where  the  geo- 
St  rophic  equili brium  given  by  Eqs.  5-r  and  5-s  do  not  hold, 
lack  of  justification  for  the  independence  of  the  trans- 

ports  from  the  zonal  circulation,  a  condition  that  is  brought 
about  by  the  dropping  of  Eqs.  5-j. 

ters  «n  Variatl°n  "  ““  "i  «*““  ^  P*rame- 

rs  can  be  determined  from  the  experimental  values  for  the  transports 

of  angular  momentum  and  sensible  hMi-  «-•  P 

laHt-nrHn  i  •  .  sensib^  heat,  respectively.  Likewise,  the 

latitudinal  variation  of  the  amplitude  of  the  stream  function  can  be 
etermined  from  the  meridional  eddy  kinetic  energy  at  Level  2,  while 
he  latitudinal  variation  of  the  amplitude  of  the  temperature  dis¬ 
turbances  may  be  deduced  from  the  experimental  measurements  of  the 
meridional  eddy  kinetic  pnonm,  r  ,  , 

eddv  kineti r  P  V  '  Measurements  of  the  zonal 

using  Eas  l  7y  my  be  Used  t0  flX  che  ««  ™">ber.  Thus, 

r,  .-s,  5-t,  and  5-u,  together  with  Eqs.  5-v  and  5-w: 


iu*  =  -  JL  )_  / da6  d6  \ 

m  afj  kB\dT"deJ 


cos  (n  +  k6)  +  ||  sin  (q  +  kfi)  J  (7 -a) 


Um  “  a  j~  kA(d06)  cos  *1  +  If  sin  q  j 


vm  ~  IF  kB  cos  C  T]  +  k 6) 


V* 

m 


_  m 


- —  kA  cos  11 
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C  7-b) 

(7-c) 

(7-d) 


f 


Using  the  definition  of  Mm,  the  identity  'u*v*'  =  (l/2){ti^9Jj|  +  ujjjv*}, 
the  angular  momentum  transport  may  be  written  as: 


o  -  1  f  da9  .  /  dg9  d6\  / _B_\2  I 
m"  4  ma2  (d9  \d0  '  d0AAf/  ) 


(7-e) 


Since 


v*2  =  ^  (v*  +  v*  +  v*  -  v*)2  =  1(^*  +  v*)2  =  -^  (SJ  +  v*)2 

4m 


the  use  of  Eqs.  7-c  and  7-d  yields: 


2  2 

r  *2  k  A 

Lvf  ]  =  — 


1  + 


8a  cos  9 


(£)]+  * 


AB 


7  7 

a  cos  9 


cos  k6  (7-f) 


The  meridional  eddy  kinetic  energy  at  Level  2  is  obtained  from  the 
geostrophic  equilibrium,  Eq.  5-s,  vdtich  yields: 


a#? 

fv!  =  m  'ST 


fv*  =  m 


90^ 

w 


or 


f(v^  +  V^)  =  m 


3(95  + 


3x 


Using  =  2tf* 


30? 

f(v*  +  v*)  =  2m  ^  =  2fv* 
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whence 


v*  +  v* 


V’ 


v*  =  -1  ;..i3  =  i  ^  =  4  jh 

2  2  2  2  m 


which,  with  Eq.  7-d,  yields: 


k2A2  =  8a2 [v*2]  cos2  0 


The  combination  of  Eqs.  7-f  and  7-g  gives  the  expression 


(*)’ 


+  2  cos  k6 


/b  \ 

[v*2] 
n  -  1 

\k?) 

evfi 

=  0 


or 


B 

Af 


J  2  ^  ' 

-  W  cos^  k6  +  — V-  -  1  - 

f  [v*2] 


cos  k6 


while  Eqs.  5-y  and  7-g  yield: 


_B_ 

Af 


2af [vi2 3  6  V  1  -  cos2  kt 


and  Eq.  7-e  with  7-g: 


doe 

de 


=  l  + 


& 


-l 


M. 


m. 


2[v^2]  cos30 


(*?  * 


G7 


Thus,  Eqs .  7-h  and  7-i  yield  B/Af  and  k6  as  a  function  of  the  statis¬ 
tics  of  the  long-term  general  circulation  involving  the  meridional 
eddy  kinetic  energies  at  Levels  1  and  2,  and  the  transport  of  sensible 
heat;  Eq.  7-j  yields  dae/de  as  a  function  of  the  transports  of  angular 
momentum  and  sensible  heat  as  well  as  the  meridional  eddy  kinetic  en¬ 
ergies  at  Levels  1  and  2.  Differentiation  of  Eq.  7-g  yields: 


1  dA  _  1 
A  cTe  7 


1  d[v*2] 

r  jt.  2  -l  d  9 

[v|  ] 


tan  0 


( 7-k) 


while  differentiation  of  the  quadratic  equation  for  (B/Af)  gives  the 
expression 


1  dB  1  dA 
B  d0  A  d0 


ctn 


0  + 


1  Af  d  [v123 

2  B  d9 


+  k  sin  k6 


d6 

d0 


_B_ 

(Af 


+  cos 


(7-«t) 


Similarly,  d6/de  can  be  obtained  either  by  differentiation  of^Eq.  7-i 

or  numerically  from  k  ^  d(k6)/d0.  Since  u?  =  (  +  um^  *  Eqs. 

4m 

7-a  and  7-b  yield  the  following  expression  for  the  zonal  eddy  kinetic 
energy  at  Level  1: 


[uf] 

[v*2]  cos20 


( 7-m) 


Where  the  term  T  is  given  by 
el 
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H 


m 


‘  2  2 

af[v|  ]  cos  0 


1  dA  1  dB 
A  31  B  d0 


1  dA  (daQ  dfi\  1  dB  dc*0 
a  ds  lie-  dej  "  b  cTe  He” 

Using  Eqs .  7-a  and  7-d,  the  expression  (l/m4)[u*v*]  is  modified  to 
the  general  form: 


m 


..  1  1  dB 

cos  ks  -  k  B  d0 


sin  k6 


( 7-n) 


4 

Using  the  expression  [u*v*]/m  =  H  /fmS  and  Eq.  7-n,  the  defini- 

3  r  mm  m 

tion  for  S  is  modified  to  yield  S(  0)  as  follows: 


S(0) 


1 

k 


/ dg6  d6  \ 

^d0  de  f 


sink6 


cos  k6  - 


1  1  dB  .  .  . 

k  B  d0  k6 


(7-o) 


Figure  17  shows  the  annual  variance  of  the  meridional  wind  com¬ 
ponent  resulting  from  both  the  transient  [v'2]  and  stationary  [v*  ] 
eddies  for  the  three  levels  as  given  by  the  most  recent  reduction  of 
available  data  (Oort  and  Rasmusson,  1971).  Figure  18  shows  similar 
data  for  the  zonal  wind  component.  Figure  19  indicates  the  behavior 
of  the  ratio  [v^]/[v*2]  as  obtained  from  the  data  in  Fig.  17,  a  ratio 
that  appears  in  Eqs.  7-h  and  7-<t. 

Figure  20  shows  the  phase  lag  parameter  as  a  function  of  the  wave 
number  in  the  range  1  sk  £  8.  The  parameter  k6  in  this  figure  is  de¬ 
rived  from  the  simultaneous  solution  of  Eqs.  7-h  and  7-i,  a  procedure 
that  requires  the  use  of  the  experimental  data  shown  in  Figs,  15,  17, 
and  19.  Note  that  for  low  latitudes  there  are  no  solutions  from  Eqs. 
7-i  and  7-h  for  the  specified  behavior  of  [v^  ]/[v^  ]  in  Fig.  19  at 
those  latitudes;  since  as  k6  -*  0  in  Fig.  20,  B/Af  becomes  very  large 
in  Eq.  7-i  but  remains  finite  in  Eq.  7-h.  This  is  a  result  that  is 
expected  from  the  breakdown  of  the  geostrophic  equilibrium,  i.e.,  Eqs. 
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c*jJw*,V  c"*1) 


NORTH  LATITUDE,  d«g 


FIGURE  17.  Experimental  Transient  Plus  Standing  Meridional  Eddy  Kinetic  Energy 
at  Levels  1,  2,  and  3  (250,  500,  and  750  mb,  respectively)  as 
Obtained  from  Oort-Rasmusson  Annual  Data  (1971) 
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NORTH  LATITUDE,  deg 


FIGURE  18.  Experimental  Transient  Plus  Standing  Eddy  Kinetic  Energy  at  Levels 
1,  2,  and  3  (250,  500,  and  750  mb,  respectively)  as  Obtained 
from  Oort-Rasmusson  Annual  Data  (1971) 
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FIGURE  20.  Value  of  k6  as  Determined  from  the  Oort-Rasmusson  Annual  Data  (1971) 
for  the  Meridional  Eddy  Kinetic  Energy  and  the  Transport  of  Sensible  Heat 
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5-r  and  5-s.  The  data  in  Fig.  20  indicate  that  d5/d0  >0  regardless 
of  latitude  and  wave  number  in  the  ranges  indicated  in  the  figure. 

The  latitudinal  variation  of  dae/de  is  indicated  in  Fig.  21,  a  result 
that  is  obtained  from  Eq.  7-j  using  the  experimental  data  in  Figs.  16 
and  17  as  well  as  the  previous  results  for  B/flf  and  d6/d0.  Since  a0  = 
J(da0/d0)d0,  a  constant  tilt  a  would  require  a  constant  value  of  da6/d6; 
thus,  Fig.  12  shows  that  the  tilt  depends  on  the  latitude.  Since 
daO/de  may  be  approximated  by  the  linear  relation: 

-0.62  0  +  0.61 

in  the  range  25  s  0  i  65°,  then 

a  =  0.61  -  0.310  +  b/0 

where  b  is  a  constant  of  integration  and  cr  as  well  as  0  are  given  in 
radians.  This  relation  then  indicates  that  the  tilt  becomes  more  pro¬ 
nounced  as  the  latitude  decreases,  a  trend  already  identified  earlier 
(Starr,  1948).  Figure  22  shows  the  latitudinal  variation  of  the  am¬ 
plitude  of  the  stream  function,  Eq.  S-v;  the  parameter  kA/a  is  obtained 
from  Eq.  7-g  and  the  experimental  data  in  Fig.  17.  Figure  22  indicates 
that  the  amplitude  of  the  stream  function  reaches  a  maximum  at  the 
middle  latitudes.  Figure  23  shows  the  latitudinal  variation  of  the 
parameter  B/Af,  which  is  derived  from  Eqs.  7-h  and  7-i  as  in  the  solu¬ 
tion  of  k6.  The  parameter  (l/A)(dA/d0)  is  shown  in  Fig.  24;  it  is  de¬ 
rived  from  Eq.  7-k  and  the  experimental  data  in  Fig.  17.  This  parame¬ 
ter  is  independent  of  wave  number  as  indicated  by  Eq.  7-k.  Figure  25 
shows  the  parameter  (l/B)(dB/d0)  as  derived  from  Eq.  7-t  and  the  ex¬ 
perimental  data  in  Fig.  19.  Using  Eq.  7-o,  the  reciprocal  of  the 
parameter  2S( 0)/a  is  as  shown  on  the  ordinate  of  Fig.  26;  this  figure 
indicates  that  S(  0)  -*  ®  at  0  “  52  deg  N,  i.e.,  at  the  latitude  where 
dH _/d0  -*  0  and  M_  ,  -•  0  as  indicated  by  Figs.  15  and  16.  Note  that 

ni  m ,  x 

these  conditions  exhibited  by  the  experimental  data  in  Figs.  15,  16, 
and  26  do  satisfy  Eq.  5-x,  which  is  one  of  the  two  fundamental  equations 
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that  allow  the  determination  of  the  large-scale  transports  of  angular 
momentum  and  sensible  heat,  the  other  fundamental  equation  being  5-q 
that  is  derived  from  the  laws  of  motion. 


FIGURE  22.  Value  of  kVa  a*  Determined  from  the  Oort-Rasmusson  Annual  Data 
(1971)  for  the  Meridional  Eddy  Kinetic  Energy  at  500  mb 
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NORTH  LATITUDE,  deg 


FIGURE  23.  Value  of  E/Af  as  Determined  from  the  Oort-Rasmusson  Annual  Datr 
(1971)  for  the  Meridional  Eddy  Kinetic  Energy  and  the  Transport  of 
Sensible  Heat 


77 


10)  as  Determined  from  the  Oort-Rasmusson 
sr  the  Meridional  Eddy  Kinetic  Energy  at 


NORTH  LATITUDE,  dtg 


0 


FIGURE  25.  Value  of  (1/B)  (dl$/d0)  as  Determined  from  the  Oort-Rasmusson 
Annual  Data  (1971)  for  the  Meridional  Eddy  Kinetic  Energy  and 
the  Transport  of  Sensible  Heat 
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FIGURE  26.  Reciprocal  of  25/ o  Parameter  as  Determined  from  the  Oort-Rasmusson 
Annual  Data  (1971)  for  the  Eddy  Zonal  and  Meridional  Kinetic 
Energy  as  Well  as  the  Transports  of  Sensible  Heat  and  Angular 
Momentum 


Figure  27  shows  the  ratio  of  the  zonal  eddy  kinetic  energy  from 
the  transient  and  eddy  motions  at  the  250-mb  level  to  the  correspond¬ 
ing  meridional  eddy  kinetic  energy  at  the  500-mb  level;  these  values 
were  obtained  from  the  most  recent  available  experimental  data  (Oort- 
Rasmusson,  1971)  and  Eq.  7-m  in  the  wave  number  range  of  interest  for 
climate.  This  figure  shows  that  a  wave  number  of  unity,  corresponding 
to  a  wave  length  equal  to  a  latitudinal  circle  with  radius  a  cos  0, 


it 
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FIGURE  27.  Ratio  of  Transient  Plus  Standing  Zonal  Eddy  Kinetic  Energy  at 
250  mb  to  the  Corresponding  Meridional  Value  at  500  mb. 
Solid  Line  Determined  from  Oort-Rasmusson  Annual  Data 
(1971).  Dashed  Lines  Determined  from  Eq.  7-m. 
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yields  the  zonal  eddy  kinetic  energy  at  the  250-mb  level.  The  dis¬ 
crepancy  at  higher  latitudes  in  the  shape  between  the  experimental  and 
analytical  curve  for  k  =  1  in  Fig.  27  may  be  attributed  to  the  shape 
exhibited  by  [u*  ]  in  Fig.  18  at  those  latitudes.  Figure  20  indicates, 
then,  that  the  phase  lag  parameter  6  61  0.1  radian  at  a  latitude  cor¬ 
responding  to  the  maximum  transport  of  sensible  heat  (Fig.  15). 


IX.  FORMULATION  OF  THE  THEORY  OF  THE  GENERAL  CIRCULATION  FOR 
CLIMATE  BASED  ON  A  VARIABLE  TILT,  PHASE  LAG,  AND 
AMPLITUDES  OF  THE  BAROCLINIC  DISTURBANCES 


The  experimental  data  for  the  transports  of  angular  momentum  and 
sensible  heat,  as  well  as  for  the  meridional  and  zonal  eddy  kinetic 
energies,  provide  the  following  results: 

1.  A  successful  test  of  the  formulation  of  the  large-scale  trans¬ 
ports,  as  given  by  Eq.  5-x,  at  the  latitude  for  maximum  trans¬ 
port  of  sensible  heat. 

2.  A  wave  number  with  a  magnitude  appropriate  for  climate  that 
includes  the  effects  of  baroclinic  disturbances  and  a  non- 
homogeneous  earth’s  surface. 

3.  The  variation  with  latitude  of  the  tilt  and  phase  lag  wave 
parameters  as  well  as  of  the  amplitudes  of  the  stream  func¬ 
tion  and  temperature  disturbances. 

These  results  indicate  that  the  assumptions  of  constant  B,  a,  and 

6  used  for  the  evaluation  of  [u*v*l/m^  in  the  constant  S  model  are  not 

mm 

verified  by  the  experimental  data.  Thus,  all  the  limitations  of  the 
constant  S  theory  described  earlier  can  be  removed  by  the  calculation 
of  a  local  S(  0)  as  given  by  Fq.  7-o,  a  procedure  that  eliminates  the 
need  for  a  rather  weak  boundary  condition  at  the  pole  to  determine  the 
constant  S.  Further,  the  available  experimental  data  can  be  used  to 
select  boundary  conditions  at  a  latitude  corresponding  to  near  the 
north  edge  of  the  Hadley  cell  (vtfiere  the  ITCZ  predominates,  i.e., 

Figs.  6,  7);  this  procedure  eliminates  the  need  for  the  use  of  a  non¬ 
existing  geostrophic  equilibrium  near  the  equator,  vrfien  considering 
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advertent  modifications  of  the  climate  of  the  earth's  surface  at  the 
middle  and  higher  latitudes.* 

The  modification  of  the  constant  S  theory  to  allow  the  determina¬ 
tion  of  the  local  S(  e)  parameter,  as  well  as  the  meridional  and  zonal 
eddy  kinetic  energies  at  Levels  1  and  2,  involves  the  use  of  expres¬ 
sions  that  relate  the  transports  of  sensible  heat  and  angular  momentum 
to  the  wave  and  amplitude  parameters  (i.e.,  Eqs.  5-y  and  7-e),  the  ap¬ 
proximation  R  M  M  , ,  and  the  eddy  kinetic  energy  at  the  250-mb 
mm  m,± 

level.  Thus,  using  a  procedure  similar  to  that  for  the  derivation  of 
Eq.  7-e,  i.e.,  the  identity  ujjjv*'  =  (l/2){u*v*  +  0*v*},  and  Eqs.  7-a 
through  7-d: 


Note  that  the  eddy  kinetic  energy  equation  could  also  be  applied  at 
Level  3;  however,  the  approximation  Mm  used  in  Eq.  7-q  yields  no 
angular  momentum  at  Level  3  as  can  formally  be  seen  from  the  following 
relations: 

[«i  =  J  ten*  -  u*Kv  -  >)]  =  -i,  [coj  -  Jjxqj  - 

4m 


VSf 

Numerical  calculations  using  tfrie  Mintz-Arakawa  model  (Gates,  iy71) 
for  two  conditions  representing  the  current  Arctic  environment  and 
’  that  of  a  postulated  Arctic  without  the  ice  indicate  no  changes  at 
these  low  latitudes  in  the  temperature,  the  geopotential  height,  and 
the  east-west  winds  in  the  lower  (i.e.,  1000- ,  800- ,  and  400-mb 
levels)  troposphere  (Warshaw  and  Rapp,  1972). 
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Expanding  the  right-hand  side  of  this  expression  and  using  the  iden¬ 
tities  for  u*v*  as  well  as  u*v*: 

mm  mm 


[U1V3]  =  T 


(7-r) 


Similarly 


Cuff}]  =  | 

or 


[(u*  +  G*)(9*  +  0*)]  =  -i*  [(Q*  +  G*)(9*  +  0*)] 

Am2  m  m/v  m  nrJ 


Cuivi]  =  r  (Hm +  V  ■  m\  =  m2«, 


m. 


(7-s) 


Thus,  the  approximation  Rm  =*■  Mm  implies  that  the  dynamics  in  the  lower 
half  of  the  atmosphere  are  relatively  unimportant  as  compared  with 
those  of  the  upper  half  of  the  atmosphere. 

The  kinetic  energy  per  unit  mass  for  the  horizontal  motion  in 
the  troposphere  is  given  by  (Lorenz,  1967) 


dK  -I 

P  1  -4  -♦  -• 

ar  =  -  o u  *  V  +  u 


(8-a) 


where  Kg  -  (1/2)  U  •  U  -  (l/2)(u^  +  v^).  Using  the  condition  of  ver¬ 
tical  hydrostatic  equilibrium  and  the  geostrophic  potential  $=  gz, 

?zp  =  «  y  =  o  u 

Thus,  Eq.  8-a  becomes 


?  J  f  U  •  F 


(8-b) 
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Equation  8-b  must  then  be  included  in  the  set  of  primitive  Eqs.  5-e 
through  5-h. 

Using  the  continuity  and  hydrostatic  equations  (5-g  and  5-h)  in 
Eq.  8-b,  making  use  of  spherical  coordinates,  taking  the  zonal  average, 
separating  the  eddy  and  circulation  terms,  and  using  a/3t  -♦  0  for  cli¬ 
mate,  Eq.  8-b  can  be  written  in  the  following  form  (Appendix  B): 


f  h  +  ^ +  +  Ip  joyw  -  cv*]  - 

+  2  ct*u.*3  -  [U*F*  +  V*F*]  -  K  iHip  +  M 

-jtap  [u*2]+f|5^  +  ^jcv]  =  0 


(8-C) 


where 


=  ^  U*  •  U*  +  [U]  •  U*  =  ^  u*2  +  ^  v*2  +  [u]  u*  +  [v]  v* 


and 


[XE]  =  \  [u*2]  +  £  [v*2’j  «K* 

When  applying  Eq.  8-c  to  the  250-mb  level,  some  of  the  terms  in  this 
equation  may  be  simplified  through  the  use  of  the  relations  for  the 
eddy  winds  (Eqs.  7-a  through  7-d),  the  geostrophic  equilibrium,  the 
boundary  conditions  for  iu  (Fig.  13),  and  the  continuity  equation 
(5-g).  Thus,  the  second  term  in  the  first  bracket  yields 

[Ke  V*]  =  £  [u*2v*]  +  i  [v*3]  +  [ux]  [u*v.i  ♦  [Vj]  [v*2] 


ee 


The  triple  correlations  in  the  right-hand  side  of  this  expression 
vanish,  as  it  can  be  seen  from  the  following: 


[ii*2v*]  =  ±  C(G*  +  u*)2(9*  +  0*)]  =  -K  [(0*  4  u*)2(9*  +  v*)] 
x  x  c  0m-5  in  m  m  m 


nr 


Expanding  the  product  (Q*  +  u£)2(G£  +  v*),  using  Eqs.  7-a  through  7-d, 
and  taking  the  zonal  average  of  each  of  the  six  individual  terms  in 
this  product,  the  results  contain  the  zonal  average  of  cubic  terms 
given  by  the  product  of  sines  and  cosines;  e.g.,  [sinc  T)  cosd  T)], 

[cosc  7]  cosd  (  71  +  k  5)  ],  [cosc  7)  sind  ( 71  +  k  6)  ],  [sinC  7)  cosd  (  7]  +  k  6)  ], 
[sin  7)  cos  71  sin  (7)  +  k6)],  etc.,  where  c  +  d  =  3  and  either  c  or  d 
take  values  0,  1,  2,  or  3.  Since  the  zonal  average  of  each  of  these 
cubic  products  vanishes,  [u*  v*]  =  0.  Similarly,  since 


[vf]  =  [(9*  +  C*3]  =  [(9J;  +  Op3] 

8m 


the  results  contain  again  four  cubic  terms  involving  the  expression 
[cosc  T|  cos^  ( T]  +  k6)].  Hence  Cv?^3  =  0.  Thus,  using  Eq.  7-s: 


CKg  v{]  =  m2  Cu,]  M  +  [Vl]  [v*2]  (8-d) 

2 

where  [v*  ]  is  given  by  Eq.  7-f.  The  third  term  in  the  first  bracket 
of  Eq.  8-c  vanishes,  as  it  can  be  seer,  from  the  use  of  the  geostrophic 
equilibrium,  i.e. 


(8-e) 


The  first  term  in  the  second  bracket  of  Eq.  8-c  may  be  written  as: 
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fp  [Kg]  W  =  fj  K*[«] 


0  "  K2CtU2]  _  2  K2[u)2] 
0  “  P2  ify 


Since  the  continuity  equation  with  the  boundary  conditions  for  u>  yields 


=  -  (am)  _  _  “b  ~  “*2  _  _ 

1  "  Wr[  "  <5  -  p2  "  p2 


one  can  write 


2  [m2]  ^  02] 

P4  =  “pT 


--ik 


[V,] 


m 


or 


3_ 

3P 


[Kg]  [tti] 


»  m  a  Cyl 

’2  a  90  m 


( 8  — f  ) 


The  condition  [F]  =  -[F]  that  was  used  in  Eqs.  5-n  and  5-o  implies 
[F1]  =  0,  hence  [F*]  =  0.  Thus,  using  Eq.  7-s  as  well  as  Eqs.  8-d 
through  8-f,  and  the  notation  3  [Xu^/dp  =  2  [X2ui|]/p4  where  X2  de¬ 
notes  Kg  2,  0*,  u*  or  v*,  Eq.  8-c  applied  at  the  upper  half  of  the 
atmosphere  yields: 


adGJm  1J  m^  m  | 


[Vg]  V*3 


+  2 


c«>j]  R 


_  jd  v*  a _ i_  ,  o 

a  *2  30  m  +  '  P4 


p  (2  j,  [u*u&]i 

T,  ‘W  -  jr  Is  \ +  2  -§rj  Cui 


-I 


tan 


[vf2] 


A  K2 3  ♦  f  le  ^  *  2  t*il  «  0  (»■») 
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Multiplying  Eq.  8-g  by  a/m,  collecting  the  two  terms  involving  M  - , 

in  |  x 

and  using  =  (1/2)  as  a  result  of  a  linear  interpolation  for  ou^, 
from  neighboring  levels,  Eq.  8-g  yields: 


h  ki  +  [vf 1} +  \  h  [ui]  m 


[vf] 


K*  ae  m  "  ^vi^  ae  m  "  ^-vi^u*  ^  sin  6 


2a 


+  HKP  <3  +  +  R  [T?<]  -  [u>*][un]  -  =  0 


mp. 


'Eg  ■  2- 


12 


2 "2 JL  1- 


2  2Jl  l1 


(8-h) 


Equation  8-h  is  the  final  form  of  the  eddy  kinetic  energy  equa¬ 
tion  at  the  250-mb  level.  The  last  bracket  contains  all  the  relevant 
vertical  transports,  Which  have  a  large  coefficient  through  the  ap¬ 
pearance  of  the  radius  of  the  earth.  A  systematic  description  of  each 
term  in  this  equation  is  given  below. 

The  meridional  circulation  at  the  250-mb  level  is  given  by  the 
2 

expression  [vn ]  =  (m  /af)  (dM  ,/dQ).  The  zonal  average  of  the  eddy 

m,i 

kinetic  energy,  x£,  can  be  defined  in  terms  of  the  wave  parameters 
and  the  amplitude  of  the  baroclinic  disturbances  with  the  aid  of 
Eqs.  7-f,  7-m,  and  7-g.  The  zonal  circulation,  [u^],  is  obtained  from 
Eq.  S-j ;  i.e.,  expanding  the  operator  d/dt  with  the  aid  of  Eq.  5-g, 
using  spherical  coordinates  for  v  •  vU,  taking  the  zonal  average, 
setting  3[v]/dt  =  0  and  using  [9]  =  0: 


m  a  plf  .  K!2L 

a  d9l  m  m  | 
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Multiplying  through  by  a/tan  0  and  rearranging  terms,  this  expression 
is  reduced  to  an  algebraic  equation  for  the  zonal  circulation  at  Level 
1 ,  i.e., 
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[ux]2  +  af  ctn  0  •  [u1  ]  +  [u*2]  -  [v^2  -  [v£2]  +  2  ctn  0  [v^ 


a  [vx] 

“5© 


+  ctn  0  -2-  j-v*^]+~  [v*u)*]ctn  0  +  ctn  0 


_ 


30 


=  0 


(8-i) 


where  3[0^]/30  is  evaluated  by  using  the  hydrostatic  equation  together 
with  a  linear  representation  of  an  emagram  giving  the  In  p  vs.  T  prop¬ 
erties  for  air;  i.e.,  using  the  relation  0  =  30p/3p  -  p30/3p  and 
Eq.  5-h: 


[0qPq]  “ 

=  -  0  .  P;  +  R  CT^  =  [«,.]  +  R  [Tl] 


Applying  Eq.  5-h  at  Level  3  with  04  =  gz4  =  0: 

[02]  =  |  R  [t3] 


hence 

[0*]  =  |  R  CT3]  +  R  [T^  (8-j ) 

2 

Since  the  constant  y  is  related  to  the  static  stability  parameter  at 
Level  2  as: 


A 


30  3 


^  ^  % 


1  2  RT  3  .  0 

’7P2"p1pln  e] 


P2 


or 


y2  =  -  |  R  p2  r2 
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where  T 2  is  the 


static  stability  parameter  at  Level  2  given  by 


[T3]  -  [TJ 
P3  “  Px 


2  tT2l 
7  P2 


then 


[T3]  -  [T^  =  j  [T2] 


(8-k) 


From  a  linear  relationship  in  an  emagram  (Kurihara,  1970): 

.  ^  ■  l1  -  -sH)  ^  +  iH  ^  <8-» 

From  Eqs.  8-k  and  8-1: 

^  2-eH# 

and 

tv-  - 1  (9  - 2  iH)  ^  - 2  i1  -  iH) [£■ 

Substituting  these  expressions  for  [T^]  and  [Tj]  in  Eq.  8-j,  and 
evaluating  9[01]/99: 


or 


9C0l3  _  flZ  10  In  2\  „  a^T23 

-rr~  ~  \r  -  n  in)  R  “5F" 


a[0x] 

se 


(8-m) 
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Note  that  the  term  3[T2 ]/ S©  in  Eq.  8-m  cannot  be  obtained  from  the  ap¬ 
plication  of  energy  Eq.  5-f  at  Level  2  to  determine  [T2 3 »  because 
[9  ]  =  0.  Since  geostrophic  Eq.  5-r  as  applied  to  the  long-term  zonal 
circulation  yields  9[  0^/90  =  -af[u  ^],  the  linear  term  in  in  Eq. 

8-i  is  approximately  equal  to  the  term  containing  9[  0^3/90  in  the  same 
equation.  Thus,  Eq.  8-i  expresses  the  degree  of  departure  from  geo¬ 
strophic  equilibrium  for  the  long-term  zonal  circulation. 

The  zonal  average1  of  the  eddy  kinetic  energy  at  Level  2,  K*,  can 
be  evaluated  with  the  aid  of  Eq.  7-g  for  [v$2]  and  the  expression  for 
[u*2].  The  latter  can  be  evaluated  from  u^2,  which  is  derived  in  a 
manner  similar  to  v*2,  i.e.,  from  geostrophic  Eq.  5-r,  u*  =  (l/2m) 
However,  Figs.  17  and  18  indicate  that  [u*2]  «  [v*2]  or  K*  «  [v*  ]. 

The  vertical  transports  involving  in  the  eddy  kinetic  energy 
Eq.  8-h  are  very  small  (Lorenz,  1967)  and  their  transient  values  ber- 
come  unmeasurable  (Oort,  1971).*  Nevertheless,  the  vertical  trans- 

A 

ports  can  be  evaluated  by  using  the  relations  uj2  =  -P2ci  =  ~5P2D  = 

A 

-4P4D,  (Appendix  A),  i.e., 


1  p4 

“  7  T 


9_  vm\ 


or  after  using  Eqs.  7-a  and  7-c, 


(u*  =  ^  ctn  0  •  cos(  T\  +  k6) 


(8-n) 


Note  that  the  large  term  a/p^  in  the  coefficient  of  the  vertical 
transports  in  Eq.  8-h  is  cancelled  out  by  its  reciprocal  in  the  co¬ 
efficient  of  Eq.  8-n.  Each  of  the  vertical  transports  in  Eq.  8-h  is 
evaluated  as  follows: 


* 

Private  communication. 
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/. 


(8-o) 


[kE2<^]  =  [u*a$][u2] 

since  the  triple  correlations  vanish  and  [9]  =  0;  the  zonal  circula¬ 
tion  at  Level  2  in  Eq.  8-o  is  given  by  the  average  of  the  values  at 
Levels  1  and  3,  i.e., 

Cu2]  =  1  ([ul]  +  [u3])  =  1  Cul]  ^  +  Xu^ j) 

where  [u^l/Cu^  <  1.  The  transport  [u*ur|]  is  evaluated  by  using 
u*  =  (1/2  m)  3^.  Hence, 


[“?•?!  =  k 


Equations  7-b,  8-n,  and  5-y  then  yield: 


[u*u)|] 


1  mp4  kHm  Ctn  9  /d«fl 
4a  af  \de 


ctn  k6  + 


1  1  dA\ 
k  A  dQ) 


(8-p) 


The  transport  [0*u)*]  is  evaluated  by  using  0*  =  (1/2)  <jS*  and  Eq.  8-n; 

the  parameter  is  obtained  from  the  geostrophic  Equation  5-s,  i.e., 
2 

=  /af  )(  a<?f/dX).  Integration  of  this  equation  with  the  aid  of 
Eq.  7-d  yields  $*  =  fA  sin  T|.  Thus: 


1 

[0*1^*]  =  -  4  —  Hm  ctn  0 


(8-q) 


The  transport  R  [T*u>* ]  or  2R  [T*uu|]  can  be  estimated  by  taking  an 
average  between  the  corresponding  values  at  Levels  0  and  2.  The 
transport  [Tgug]  vanishes  because  ug  =  0,  while  the  transport  [Tgug] 
is  evaluated  by  using  0*  =  RT*  and  Eq.  8-n.  The  result  is: 
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R  =  [tf* u>*]  =  0 


(8-r) 


Since  [T*u>*]  =  [T*u$]  =  0,  its  average  [T*uif]  =0.  The  transport 
[v*u^]  is  evaluated  by  using  v*  =  (^)  v£  and  Eq.  8-n;  the  result  is 
then  as  follows : 

[vjtuj]  =  •  Hm  ctn  0  .  ctn  k6  (8-s) 

The  theoretical  determination  of  the  meridional  and  zonal  kinetic 

energy  terms  at  Levels  1  and  2  in  Figs.  17  and  18,  as  well  as  the 

parameters  in  Figs.  20  through  27,  can  be  obtained  from  the  numerical 

integration  of  the  foregoing  system  of  equations.  These  numerical 

integrations  can  start  with  initial  conditions  at  low  latitudes  that 

are  determined  from  these  figures.  The  initial  conditions  involve, 

for  example,  the  terms  that  appear  in  Eq.  7-q,  i.e.,  Mm  Hm,  dae/de, 

d6/d0,  6,  dlnA/d9,  and  dlnB/dfl.  The  basic  Equation  5-q  yields  then 

dMm  ^^/de  for  a  given  value  of  the  eigenvalue  parameter  y  ;  likewise, 

Eq.  5-x  yields  df^/dQ  from  the  initial  values  for  Hm,  Mm  1,  and  S. 

The  problem  is  then  reduced  to  the  determination  of  four  derivatives, 

i.e.,  d2of0/d02,  d26/d02,  d2lnA/d02,  and  d2lnB/d0,  which  are  obtained 

from  the  simultaneous  solution  of  four  ordinary  differential  equations. 

These  equations  are  given  by  Eq,  8-h  and  by  three  equations  obtained 

through  differentiating  Eqs.  7-e  and  7-q  once  and  differentiating  Eq. 

5-y  twice;  the  second  derivative  of  Hm  which  then  appears  can  be 

evaluated  from  the  differentiation  of  Eq.  5-x  or  5-z.  Similarly,  the 

second  derivative  of  Mm  which  appears  in  d[v^]/d0  in  Eq.  8-h,  can 

be  obtained  from  Eq.  5-q.  The  numerical  integrations  can  be  performed 

2 

until  the  eigenvalue  parameter  y  is  found  to  satisfy  the  boundary 
condition  -*  0  as  0  -•  tt/2. 

The  solution  of  the  foregoing  system  of  equations  yields  the  dy¬ 
namics  induced  by  the  heating  rates  of  Eq.  1;  note  that  they  are  de¬ 
coupled  from  the  transport  of  water  vapor  by  two-dimensional  turbu¬ 
lence.  The  converse  is  not  true,  however,  for  the  large-scale  transport 
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of  water  vapor  is  coupled  to  the  dynamics  through  the  continuity  equa¬ 
tion  for  water  vapor.  Equation  5-c  applies  to  either  a  gas  or  a 
liquid;  however,  it  is  usually  assumed  that  liquid  water  falls  out 
im  lediately  upon  forming  from  condensation.  Thus,  the  continuity 
equation  for  water  vapor  is  given  by  a  modification  of  Eq.  5-c  to  in¬ 
clude  the  effects  of  condensation  and  diffusion  from  mainly  the  ver¬ 
tical  gradient  of  moisture.  Since  differentiation  of  Eq.  5-c  yields 
the  expression  dp/dt  =  -p  7  •  V  for  air,  the  respective  equation  for 
water  vapor  is  given  by  (e.g.,  Washington  and  Ka Sahara ,  1970): 

^  =  -Pw  V  .  V  +  Mw  +  pEw  (9-a) 

where  pw  denotes  the  density  of  water  vapor,  Mw  the  rate  of  condensa¬ 
tion  of  water  vapor  per  unit  volume,  and  Ew  the  rate  of  change  of 
water  vapor  content  per  unit  mass  due  to  mainly  the  vertical  diffu¬ 
sion  of  water  vapor.  Equation  9-a  can  be  written  in  terms  of  either 
the  specific  humidity  or  the  mixing  ratio,  which  are  defined  as  the 
amounts  of  water  vapor  per  unit  mass  of  moist  or  dry  air,  respectively. 
Since  the  specific  humidity  differs  only  slightly  from  the  mixing 
ratio,  the  use  of  either  one  is  approximately  interchangeable.  If 
the  specific  humidity  is  denoted  by  q,  i.e.  (e.g.,  Willett  and  Sanders, 
1959): 


q 


0.622  e 
p  -  0.378e 


m  0.622  - 
P 


(9-b) 


where  e  denotes  the  partial  pressure  of  water  vapor,  Eq.  9-a  and  the 
use  of  the  relation  q  dp/dt  =  -pq  V  •  V  yields 


w 


(9-c) 


which  in  the  coordinate  system  of  Eq.  5-g  (e.g.,  Holloway  and  Manabe, 
1971)  yields: 
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(9-d) 


m  9  [vq]  3 
a  "58  m  "Sp 


[iuq]  = 


+  ^ 


Applying  Eq.  9-d  to  the  atmospheric  column  in  Fig.  1,  and  denoting  the 
vertical  average  for  each  term  by  X2  =  jJ  X  d(p/p4)  where  X  denotes 
[vq],  [wq],  [Mw/p],  and  [Ew]: 


m  a  0^2 ] 

a  39  m 


[Ew  ^ 
2 


(9-e) 


Note  that  the  first  term  in  the  right-hand  side  of  Eq.  9-e  is  related 
to  the  heat  released  from  condensation,  C,  in  Eq.  1  as  follows: 


or 


C  = 


(  9-f  ) 


vdiere  L  denotes  the  latent  heat  of  condensation 
heating  rate  per  unit  mass,  i.e.,  =  L  Mw/p, 


.  Also,  if  denotes  the 
Eq.  9-f  can  be  written: 


(9-g) 


Similarly,  the  rate  of  evaporation  per  unit  mass  over  the  oceans  that 
is  injected  into  the  atmospheric  column  can  be  related  to  [E  _].  Note 

w,  4 

that  the  transport  of  specific  humidity  as  given  by  Eq.  9-e  can  be  ex¬ 
pressed  as  the  difference  of  the  rates  of  evaporation  from  the  earth 
and  precipitation  upon  the  earth  (i.e.,  Fig.  6);  this  result  is,  of 


96 


course,  in  agreement  with  that  derived  by  using  Eq.  6  with  X  =  q 
(Lorenz,  1967).  The  transport  of  water  vapor  itself  in  Eq.  9-e  can 
be  evaluated  directly  as  a  function  of  latitude,  since  the  vapor  pres¬ 
sure  is  a  function  of  temperature  (e.g.,  Petterssen,  1956);  thus,  from 
Eq.  9-b,  q*  becomes  proportional  to  .  It  is  then  seen  that  through 
the  transport  of  water  vapor,  the  dynamics  of  the  troposphere  control 
the  amounts  of  evaporation  from  the  earth  and  precipitation  upon  the 
earth . 

The  application  of  the  foregoing  equations  to  the  problem  of  ad¬ 
vertent  modification  of  the  climate  of  the  earth’s  surface  at  middle 
and  high  latitudes  involves  the  use  of  a  procedure  for  iterating  T2( e) 
for  a  given  advertent  modification  of  the  earth’s  surface  albedo  at 
those  latitudes.  The  occurrence  of  evaporation  and  condensation  in 
the  atmosphere  is  caused  primarily  by  the  changing  capacity  of  the  air 
for  water  vapor,  which  depends  entirely  upon  the  air  temperature  (e.g., 
Willett  and  Sanders,  1959).  Thus,  the  use  of  an  arbitrary  initial 
T2(6),  together  with  Eqs.  1  and  2,  generates  the  function  qf ( 6) ,  which 
in  turn  yields  the  dynamics  as  well  as  the  humidity  transports.  The 
subsequent  iteration  of  T2( 0)  may  then  proceed  until  Eq.  9-e  is  satis¬ 
fied.  Thus,  this  iteration  procedure  starts  with  the  use  of  an  a  priori 
modification  of  the  parameterization  of  the  troposphere  heating,  say, 
for  example,  at  the  middle  and  high  latitudes  as  the  result  of  the  ad¬ 
vertent  modification  or  the  Arctic  Ocean's  surface  albedo  (Fig.  14). 
Therefore,  this  a  priori  modification  must  account  for  the -effects 
present  in  Eqs..  1  and  2,  i.e.,  the  latent  heat  C,  the  local  albedo  of 

the  atmosphere  (S  )  and  the  earth’s  surface  (S*),  etc. 

61 

It  is  emphasized  that  the  main  assumption  in  the  foregoing  system 
of  equations  based  on  the  current  climate  is  the  important  role  of  the 
dynamics  of  the  upper  half  relative  to  those  of  the  lower  half  of  the 
troposphere,  i.e.,  M  “■  M  or  M^/M^  «  1.0  (Table  3).  Because  the 

transport  of  angular  momentum  in  the  upper  half  of  the  troposphere 
controls  the  meridional  cell  structure  of  the  general  circulation, 
which  prevails  even  during  the  two  extremes  of  the  interglacial  and 
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glacial  cycles  of  climate  during  the  current  Pleistocene  Epoch,  the 
assumption  concerning  the  relative  role  of  the  dynamics  in  the  upper 
half  of  the  troposphere  becomes  applicable  to  the  general  circulation 
of  the  atmospheric  air  for  any  climate. 
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X.  RESULTS  AND  CONCLUDING  REMARKS 

Considerations  of  the  feasibility  of  advertent  climate  modifica¬ 
tion  at  the  higher  latitudes  of  the  United  States  involve  the  answers 
to  four  key  subquestions,  which  are  identified  in  Section  II.  This 
paper  seeks  to  answer  the  first  or  initial  subquestion  through  the 
formulation  of  a  model  for  the  general  circulation  for  climate  that 
emphasizes  physical  phenomenology.  This  model  is  based  on  a  variable 
(i.e.,  as  a  function  of  latitude)  tilt  from  the  meridian  for  both 
the  flow  and  temperature  disturbances,  a  variable  phase  lag  between 
the  flow  and  temperature  disturbances,  and  a  variable  amplitude  of 
these  baroclinic  disturbances.  Results  from  this  formulation  are  as 
follows : 

1.  The  use  of  this  model  utilizing  experimental  data  for  the 
long-term  statistics  of  the  general  circulation  (i.e.,  the 
annual  zonal  averages  of  the  transports  of  angular  momentum 
and  sensible  heat  as  well  as  the  meridional  and  zonal  eddy 
kinetic  energies  at  the  250-  and  500 -mb  levels)  allows  the 
semi-experimental  determination  of  parameters  that  are  dif¬ 
ficult  to  measure  directly,  i.e.,  the  wave  and  amplitude 
parameters  of  the  baroclinic  disturbances,  the  wave  number 
appropriate  for  climate,  and  the  vertical  transports  (i.e., 
Figs.  20  through  27,  and  Eqs .  8-o  through  8-s).  These  re¬ 
sults,  for  example,  indicate  quantitatively  that  the  tilt  of 
the  baroclinic  disturbances  becomes  predominant  as  the  lati¬ 
tude  decreases,  which  is  in  agreement  with  previous  qualita¬ 
tive  experimental  observations  (Starr,  1948).  Also,  the 
wave  number  appropriate  for  climate,  and  including  the  effects 
of  baroclinic  disturbances  as  well  as  of  a  nonhomogene ous 
earth’s  surface,  is  of  the  order  of  unity. 
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2.  There  is  good  agreement  between  the  experimental  value  of  the 
middle  latitude,  at  which  the  large-scale  transport  of  angular 
momentum  vanishes  as  the  transport  of  sensible  heat  retches 

a  maximum  value  (Pigs.  15,  16),  with  the  corresponding  lati» 
tude  value  that  is  determined  from  the  semi-experimental 
values  of  the  wave  and  amplitude  parameters  of  the  baroclinlc 
disturbances  (Fij.  26,  Eqs.  7-o  and  5-x).  This  agreement  con¬ 
stitutes  then  a  successful  local  test  of  the  basic  formulation 
of  the  baroclinic  disturbances  for  the  determination  of  the 
large-scale  transports  of  angular  momentum  and  sensible  heat 
by  two-dimensional  turbulenc?  (i.e.,  Lqs.  5-v  and  5-w) . 

3.  The  model  bast  I  on  variable  wave  and  amplitude  parameters  of 
the  baroclinic  disturbances  indicates  that  the  large-scale 
transports  of  angular  momentum  and  ser^ible  heat  by  two- 
dimensional  turbulence  do  depend  on  the  zonal  circulation 
(i.e.,  Eqs.  8-h  and  8-i),  a  result  that  resolves  the  puzzle 
of  the  constant  S  theory  which,  by  implicit  assumption, 
yields  the  result  that  the  transports  are  independent  of 
the  zonal  circulation  (Smagorinsky,  1964).  The  model  based 
on  variable  wave  and  amplitude  parameters  of  the  baroclinic 
disturbances  indicates  further  that  the  transports  of  angular 
momentum  and  sensible  heat  are  also  coupled  to  the  vertical 
transports  as  given  by  Eqs.  8-h  and  8-o  through  8-s. 

4.  The  model  based  on  variable  wave  and  amplitude  parameters 
of  the  baroclinic  disturb  ncoa  requires  only  one  eigenvalue 
parameter,  which  Involves  the  static  stability  of  a  tropo¬ 
spheric  column;  this  result  implies  that  there  is  no  need, 
as  in  the  constant  S  model,  for  the  use  of  a  rather  weak 
boundary  condition  at  the  north  pole  in  order  to  establish 
the  magnitude  of  the  parameter  S  through  most  of  the  lati¬ 
tudinal  range  of  the  northern  hemisphere. 

The  model  based  on  variable  wave  and  amplitude  parameters  of 
the  baroclinic  disturbances  may  also  modify  an  important  conclusion 
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that  is  derived  from  the  constant  S  model,  which  indicates  that  the 
same  multiple  of  the  transports  of  angular  momentum  and  sensible  heat 
ar  well  as  of  the  poleward  heat  flux  satisfies  the  fundamental  oqua- 
fions  of  the  model  (i.e.,  Eqs .  5-q  and  5-x  with  S  a  constant).  This 
would  then  moan  that  the  transports  scale  directly  with  the  magnitude 
of  the  poleward  heat  flux,  i.e.,  that  only  the  shape  of  the  poleward 
heat  flux  function  (Fig.  la)  determines  the  two  or  three  meridional 
coll  structures  (e.g.,  Fig.  5)  corresponding,  respectively,  to  the  two 
extremes  of  interglacial  or  glacial  climates  of  middle  and  hir,h  lati¬ 
tudes  during  the  current  Pleistocene  Epoch.  Thus,  the  shape  of  the 
poleward  heat  flux  function  would  determine  the  type  of  the  climate, 
while  the  amplitude  of  the  poleward  heat  flux  would  fix  the  intensity 
of  the  particular  type  of  climate.  However,  the  model  based  on  vari¬ 
able  wave  and  amplitude  parameters  indicate*  that  S  is  a  function  of 
latitude,  which  depends  on  the  meridiont*  and  zonal  eddy  kinetic  energy 
as  well  as  on  the  vertical  transports;  hence,  the  foregoing  conclusion 
from  the  constant  S  model  may  become  invalid. 

The  formulation  of  the  model  based  on  variable  wave  and  ampli¬ 
tude  parameters  of  the  barocllnlc  disturbances  remains  to  be  tested 
for  the  current  climate.  This  test  requires  the  following  two  steps: 

O'  to  obtain  for  latitudes  greater  than  about  30  degrees  (Fig.  20) 
the  theoretical  values  of  the  transports  of  sensible  beat  and  angular 
momentum  (Figs.  15,  16)  based  on  S(0)  as  well  as  the  theoretical  values 
of  the  meridional  and  zonal  eddy  kinetic  energy  at  the  250-  and  500-mb 
levels  (Figs.  17,  18);  and  (2)  to  include  (at  the  same  latitudes)  com¬ 
parisons  of  theoretical  values  for  the  statistics  of  the  water  vapor 
(Eq.  9-e)  with  corresponding  experimental  values  from  the  Oort-Rasmusson 
data.  This  latter  consideration  might  require  a  thermodynamic  param¬ 
eterization  of  the  moisture  flux  based,  for  example,  on  seasonal  data 
for  the  current  climate. 

The  model  based  on  variable  wave  and  amplitude  parameters  of  the 
baroclinic  disturbances  could  then  be  applied  to  conditions  of  an 
Arctic  Ocean  without  sea  ice,  which  corresponds  to  the  Interglacial 
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environment  of  the  current  Pleistocene  epoch.  This  procedure  would 
use  known  boundary  conditions  at  low  latitudes  and  properly  modified 
functions  of  the  poleward  heat  flux  at  middle  and  high  northern  lati¬ 
tudes  (Pig.  14);  such  calculations  would  thus  provide  a  test  of  the 
iteration  procedure  for  the  solution  of  the  problem  of  advertent  modi¬ 
fication  of  the  earth's  climate  at  middle  and  high  latitudes. 

The  model  based  on  variable  wave  and  amplitude  parameters  of  the 
baroclinic  disturbances  could  next  be  used  to  answer  the  second  and 
th^nd  subquestions  identified  earlier  in  the  key  questions  in  ad¬ 
vertent  climate  modification  (Section  II),  i.e.,  the  predictability 
of  changes  of  the  climate  of  the  earth’s  surface  at  middle  and  high 
latitudes  of  the  northern  hemisphere,  and  the  sensitivity  of  the  mag¬ 
nitude  of  the  long-term  statistics  of  the  general  circulation  to 
advertent  modifications  of  parameters  such  as  the  earth’s  surface 
albedo,  etc.  The  transitive  or  intransitive  nature  of  the  advertent 
changes  of  climate  at  middle  and  high  northern  latitudes  could  be 
investigated  numerically  with  this  model  by  attempting  to  find  an 
extra  solution  for  the  conditions  of  an  Arctic  Ocean  without  sea  ice. 
The  sensitivity  of  the  long-term  statistics  of  the  general  circulation 
could  subsequently  be  investigated  numerically  with  this  model  by 
using  parametric  perturbations  of  the  poleward  heat  flux  at  the  middle 
and  high  northern  latitudes. 

Preliminary  investigations  of  the  fourth  subquestion  concerning 
the  technology  for  advertent  climate  modification  at  the  higher 
northern  latitudes  are  given  in  Section  VI.  In-depth  investigations 
of  this  subquestion  must  await  positive  answers  to  the  second  and 
third  subquestions. 
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APPENDIX  A 


DERIVATION  OF  ZONAL  MOMENTUM  AND  ENERGY  EQUATIONS 
AND  S  WAVE -PARAMETER 

This  appendix  gives  a  more  detailed  but  unified  presentation  of 
all  the  assumptions  that  are  used  in  (1)  the  derivation  of  the  two 
zonal  Equations  5-k  and  5-1  as  well  as  the  energy  Equation  5-m  from 
the  set  of  primitive  Equations  5-e  through  5-h,  and  (2)  the  use  of  the 
mechanisms  of  the  tilt  and  phase  lag  of  the  baroclinic  disturbances 
(Fig.  12,  Section  VII)  which  leads  to  Eq.  5-x. 

1.  ZONAL  MOMENTUM  AND  ENERGY  EQUATIONS 

The  scalar  Eq.  5-i  can  be  written  in  the  H creator  (x,  y)  coor¬ 
dinates  instead  of  the  spherical  (X,  9)  coordinates  through  the  trans¬ 
formations  dx  =  a  d\  and  dy  =  a  d 0/cos  t.  The  zonal  and  meridional 

winds  in  the  Mercator  plane  then  become  u  =  dx/dt  =  a  dx/dt  and  v  = 

.  m  m 

dy/dt  or  vm  -  a  cos  8  d9/dt,  while  the  relationship  between  the 

Mercator  winds  and  actual  winds  becomes  u  =  u  cos  0  and  v  =  v  cos  0. 

Substituting  this  expression  for 'u  in  Eq.  5-i,  the  Mercator  zonal  wind 

is  given  by  (Sn egorinsky,  1958): 


u  v 
m  m 


sin  8  +  fv 

m 


J-  |  i 

cos2e  ax  cos  6 


(A-l) 


where  the  operator  d/dt  is  given  by  3/3t  +  umd/3x  +  v^/Sy  +  u>d/dp. 
After  expanding  the  operator  d/dt  in  Eq.  A-l,  using  the  notation 
m  =  cos  0,  the  transformation  dy  -  am  d@,  and  noting  the  identity 


m  dy 
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UmVm 

m  m 


sin  8  = 


Preceding  page  blank 
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the  local  derivative  of  the  zbnal  wind  becomes 


du 


-u 


m 


m  dx 


sin  0  -  m 


2 

Sx 


+  nur 


(A-2) 


The  application  of  Eq.  A-2  at  Levels  1  and  3  (Fig.  13,  Section  VII) 
requires  the  evaluation  of  the  term  uidu^/dp  at  these  levels;  with  the 
notation  X  =  X1  -  X^,  the  term  at  Level  1  becomes 


A 


U)-,U 
l  m 

Ap 


The  value  of  is  obtained  from  the  continuity  Eq.  5-g  after  its 
transformation  to  the  Mercator  plane;  denoting  the  horizontal  diver¬ 
gence  V  •  U  by  D,  using  spherical  coordinates  in  v  •  U  and  the  fore¬ 
going  transformations  for  both  the  coordinates  and  the  winds,  one 
obtains : 


D  =  v 


ff  =  IS  +  J. 
a  ae  m / 


Oil) 

dp 


hence 


Similarly 


D 


3 
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with  the  notation  X  =  +  X^,  B  =  0,  while  D  =  2D^.  Interpolating 

linearly  for  and  ou^  from  neighboring  levels,  =  u)2/2.  Thus 


U) 


1 


A 


Using  the  identity 


u  =  2  u 
m 


m. 


“  u 


m 


Substituting  this  result  in  Eq.  A-2  at  Level  1,  and  noting  that 

u  v  sin  0/a  =  -mu  v  dm”1/dy,  Eq.  A-2  at  Level  J  can  be  put  in 

m^ 

the  form 


Similarly,  Eq.  A-2  when  applied  to  Level  3  yields 


Equations  A-3  and  A-4  are  similar,  except  for  the  change  in  sign  of 
Du  /4,  which  comes  about  from  the  use  of  u  instead  of  u  ,  i.e., 

m  m3  IU1 

the  identify  -u  =  2u  -  u  .  Adding  Eqs.  A-3  and  A-4  one  obtains 
m  m 


while  subtracting  Eq.  A-4  from  Eq.  A-3  yields 


at 


-S'- ■‘■i#)* 


Taking  the  zonal  average  of  these  equations,  defined  by 


[X] 


2rra 

■*/ 


X  dx 


one  gets 


at 


4  a  ^UmVm]  .  „rp 

u  X  J 
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^  +  ^  ♦  f[0_j  +  m[F„] 
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J  m 


(A-5) 


(A-6) 


where  in  the  derivation  of  Eq.  A-5,  the  term  [^ml  drops  out,  because 
of  D  =  0  as  well  as  the  boundary  condition  vm  =  0  at  the  equator  and 
the  poles  at  both  levels  and  at  any  time,  t.  By  using  the  identity 
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xT!  =  i  (x  9  +  x  v  ) 

mm  2  '  m  m  m  m 


m  mm 


and  the  definitions  X  =  [X  ]  +  X  =  [X  ]  +  X*  one  gets  the  fol- 

m  mJ  m’  m  m  m  3 

lowing  result 


[uV]  =  i  [u  ][v  ]  +  [u*v*] 
m  mJ  2  mJ  mJ  m  mJ 


A  _  A 


Hence,  Eq.  A-5  becomes  (Smagorinsky,  1964) 


it 


4 

m  +  m[F  ] 
ay  L  x J 


(A-7) 


(Eq.  5-k,  Section  VII) 


where  Mm  denotes  a  quantity  proportional  to  the  nonlinear  eddy  flux  of 
angular  momentum,  i.e., 


M 


m 


[U*v*]/m4 


Similarly,  by  using  the  identity 


(X  Q 
m  n 


+  X  v  ) 
m  m 


one  gets 


[u^v-,  =  i  [G  ][v  ]  +  [u*v*] 
L  m  m-  2  L  mJL  mJ  LmmJ 


and  Eq.  A-6  becomes 


It - T  8y  T  m  ay  —TT  +  —  +  f['-J  + 


m[Fx] 
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From  the  condition  B  =  0,  and  the  definitions  D  =  [D ]  +  D*, 

u  =  [u  ]  +  u*  ,  the  third  term  in  the  right-hand  side  of  this  equa- 
m  u  m.  ’ 

tion  becomes 


Expanding  the  first  term  in  the  right-hand  side  of  d[um]/dt  and  using 
the  definition  D  =  dum/dx  +  d(vm/m2)/'ay  together  with  [9m]  =  0  in 
[D][um],  the  expression  for  d[um]/dt  reduces  to  the  following  form: 


d[Gj 

m 

dt 


T  «m]  £ 


[OL] 

m 


-  m 


m 


4  a 

"Sy 


m  mJ  +  [D*u* ]  +  f[vm]  ♦  m[Fx] 


m 


or  rearranging  terms 


mJ  _ 


dt 


m 


.  ,  2  »,  tD*uJ] 

+  »  l-«  jT-  +  - J— 

7  in 


+  m[FxJ 


Using  the  condition  [9  ]  =  0,  and  the  definition  of  the  vertical  com- 

m  2  2 

ponent  of  the  relative  vorticity  C  =  ^/dx  -  m  )/9/  which 

yields  [(;]  =  -m2  d(  [um]/m2)/dy ,  the  expression  for  d[um]/dt  may  be 

written  as  follows  (Smagorinsky,  1964) 


d[u  ]  ^  \  2 

r-  =  WJ[f  +  cl  +  "> 


dt 


„  dM  [D*u*  ] 
2  m  .  L  mJ 

"  m  W  +_ 2~ 


m 


+ 


(A-8) 


(Eq.  5-4»  Section  VII) 


Equations  A-7  (5-k)  and  A-8  (5-0  have  been  derived  from  the  zonal 
component  of  Eq.  5-e  (Section  VII)  and  the  continuity  Eq.  5-g  (Section 
VII)  as  applied  to  Levels  1  and  3.  The  energy  Equation  5-f  (Section 
VII)  together  with  the  hydrostatic  relation  5-h  (Section  VII)  are  ap¬ 
plied  at  Level  2  (i.e.,  Fig.  1,  Section  III).  Since  Eq.  5-f  can  be 
rewritten  as 
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d  In  T  _  R  u» 
-3F--  cpp 


and  Che  definitions  of  w  and  Che  poCential  temperature,  0  ,  where 
0p  =  (P4/P^R  P  *  T>  y^ld 


d  In  T  _  R  uj  d  ^  6i 
cfii  Cp  p  +  dc 


the  energy  Equation  5-f  at  Level  2  becomes 


d 

3t 


In  9. 


Since  Eq.  S-h  (Section  VII)  yields  p2  =  -p/5  =  p2/5  or  T2  =  5/R,  the 

above  expression  becomes  with  the  use  of  the  definition  of  6  : 

P 


Note  that  this  expression  with  5  =  RT2  is  compatible  with  Eq.  5-f  at 
Level  2  only  because  the  term  involving  the  "vertical"  velocity 
(i.e.,  uj  =  dp/dt  =  -pgw)  is  small  compared  with  Q<,/c  .  Using  the 

^  ^  HP  Q 

operator  d/Ot  in  the  Mercator  plane  and  2U0  =  +  Uj  =  U,  the  energy 

equation  at  Lavel  2  becomes 


A 

Using  the  condition  B  =  3um/ax  +  m  ^n/m2^  =  °>  the  foregoing 
equation  can  be  put  in  the  following  form: 


d$  _ 

at 
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The  third  term  in  the  right-hand  side  may  be  evaluated  with  the  aid  of 
the  continuity  equation  and  the  definition  of  9p}  this  latter  condi¬ 
tion  yields  3<ln  6  2)/&p  x  Thus*  w2  n'*y  **  CX‘ 

pressed  as  where  =  ( 1/2 energy 
equation  at  Level  2  becomes 


Taking  the  zonal  average,  using  [0n]  *  0,  and  the  definition  of  D  to 
gether  with  the  fact  that  y2  is  constant,  the  foregoing  equation  be¬ 
comes  (Smagorinsky,  1964): 


to*1  (A'9) 

(Eq.  5-m,  Section  VII) 

where  H  denotes  a  quantity  proportional  to  the  nonlinear  eddy  flux  of 
m 

sensible  heat,  i.e., 


PS  9*3 


2.  THE  S  WAVE  -PARAMETER 

A  feature  of  the  circulation  in  middle  and  higher  latitudes  that 
is  almost  as  prominent  as  hydrostatic  equilibrium  is  geostrophic 
equilibrium,  i.e.,  the  approximate  balance  between  the  Coriolis  force 
and  the  horizontal  pressure  gradient  force  in  Eq.  5-e  (Section  VII). 
The  geostrophic  equilibrium  then  yields  -f£  x  U  =  or  equivalently, 
U  *  (l/f)£  x  Using  spherical  coordinates  for  70,  the  relationship 

between  the  Mercator  winds  and  actual  winds,  and  separating  the  cir¬ 
culation  and  eddy  motions,  one  obtains: 
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( A-10) 


O  a 

m  djr  _  m  bf 

T  "5v - 17  'ST 


(Eq.  5-r,  Section  VII) 


m 


2 

m  djr 

T  *Sx 


=  a7  "Sx 


(A-ll) 


(Lq.  5-s,  Section  VII) 


The  differentiation  of  the  expression  for  Hm»  use  of  the  condition 
D*  =  0  in  the  product  and  use  of  Eqs.  A-10  (5-r)  and  A-ll 
(5-s)  yields 


1  .  tuft]  tftfti 


T  "5y~ 


2m 


2m 


Expanding  M  =  [u*v*]/m  ,  the  foregoing  expression  becomes 
*  m  mm 


(A-12) 


The  condition  B*  =  0  yields  the  stream  function  ♦*  in  the  following 
form: 


ft  -  -1"2  #  ■  - 5  IF  <A‘13) 

(Eq.  5-t,  Section  VII) 

ft’-^-^lf  (A'14) 

(Eq.  5-u,  Section  VII) 
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Thus,  the  two  spectral  modes  (U*m  and  Q£)  of  the  eddy  winds  are  given 
by  Eqs.  A-10,  A-ll,  A-13,  A-14  and  the  analytical  description 
of  the  baroclinic  disturbances,  i.e., 


♦*  =  A  sin  T)  (A-15) 

(Eq.  5-v,  Section  VII) 

$*  =  B  sin  (71  +  k6)  (A-16) 

(Eq.  5-w,  Section  VII) 

where  71  =  k(\  -  ere).  Substituting  these  results  in  the  expression 

for  H one  gets; 
m 


Hm  =  [sin  (71  +  kfi)  cos  71]  =  sin  k6  (A-17) 

(Eq.  5-y,  Section  VII) 

Assuming  that  B,  or,  and  6  in  Eq.  A-16  (5-w)  are  nearly  constant, 


[Uy  ■  =  [cos  (71  +  kfi)  cos  71]  = 

rn  a  rm 


2a“fm 


( A-18) 


Note  that  the  conditions  of  constant  B  and  ft  used  in  the  derivation 
of  Eq.  A-18  imply  from  Eq.  A-17  that  ^  is  proportional  to  A.  Com¬ 
bining  Eqs.  A-17  and  A-18: 


Hm 
fmS 


m 


where 


.  a  tan  kfi 

T3H 


(A-19) 


(A-20) 
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Combining  Eqs.  A-12  and  A-19: 


( A-21) 

(Eq.  5-x,  Section  VII) 
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APPENDIX  B 


DERIVATION  OF  EDDY  KINETIC  ENERGY  EQUATION 

This  appendix  gives  a  more  detailed  derivation  of  the  eddy  kinetic 
energy  Equation  8-c  (Section  IX)  fiom  Eq.  8-b  (Section  IX),  i.e., 


5 

dt 


-U  •  7  0  +  U  •  F 


(B-l) 

(Eq.  8-b,  Section  IX) 


which  must  be  added  to  the  set  of  primitive  Equations  5-e  through  5-h 
(Section  VII).  Multiplying  the  continuity  Equation  5-g  through  by  -0 
and  adding  the  result  to  B-l  (8-b): 


dKe 

at" 


=  -V 


<#u  .  „  |2  +  cr .  ? 


Using  Equation  5-h,  to  eliminate  90/ Bp  yields 


— •  90ou  RTti)  r* 

ar  =  -» •  -  f — r +  u ' F 


(B-2) 


Since  Of  =  [U]  +  U*,  the  kinetic  energy  Kg  may  be  expressed  as 


follows : 


Xe  =  Xz  +  XE  (B-3) 

vrtiere  Kz  and  K^,  denote,  respectively,  the  zonal  and  eddy  kinetic  en¬ 
ergies,  which  are  defined  as  follows: 


Preceding  page  Mink 
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Kz  =  |  [ff]  •  [ff]  =  \  [u]2  +  |  [v]2 


and 

Ke  =  £  t?v  •  U*  +  [U]  •  U*'  =  2  u*2  +  J  v*2  +  [u]  u*  +  [v]  v* 

Substituting  Eq.  B-3  in  B-2,  the  individual  time  rate  of  change  of  the 
eddy  kinetic  energy  is  given  by 


dKz  dKE 

3F"  +  Ht 


-v 


( B-4) 


where  dK^/dt  =  [U]  •  d  [U]/dt.  Multiplying  the  continuity  Equation 
5-g  through  by  X,  adding  the  result  to  the  expression  for  dX/dt, 
letting  X  =  KE>  using  spherical  coordinates  for  v  •  XU*,  and  taking 
the  zonal  average: 


[KPv] 


d  [v  1  -  d  r-y  -1  III  ^  E 

at  LKE]  -  -St  [KE]  +  a  -^9 


m 


+  -Sp  [ke^ 


which,  with  the  zonal  average  of  the  right-hand  side  of  Eq.  B-4, 
yields : 


dK 


Z  ,  d  r„  .  m  3  ^(KE  + 


at^  +  ^  CKE]  +  ? 


m 


+  ^  [(KE  + 


=  -  “  [Tu>]  +  [U  •  F] 


(B-5) 


Separating  the  eddy  and  circulation  terms  in  Eq.  B-5,  i.e.,  using 
X  =  [X]  +  X*  where  X  denotes  successively  the  geopotential,  the  merid¬ 
ional  wind,  the  individual  time  ra’.e  of  change  of  the  pressure,  and 
the  temperature: 
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_a_ 

at 


[ke]  + 


m  _9_ 

a  90  1  m 


CKev*] 


"  I 


+  {[KE][a»]  +  [KEa J*]  +  +  T  =  0 


(B-6) 


where  the  term  T  denotes 
e2 

Te2  =  (dKz/dt)*  +  f  "  [u*F'x  +  v*Fe] 

and  (dX^/dt)*  the  eddy  terms  in  dX^/dt  or  [U]  •  d  [U]/dt,  i.e.,  in 

[u]  d  [u]/dt  +  [v]  d  [v]/dt.  The  expression  for  (dX^/dt)*  can  then 
be  obtained  by  taking  the  zonal  average  of  Eqs.  5-i  and  5-j  (Section 
VII),  expanding  the  individual  time  derivatives  in  a  manner  similar 
to  that  of  d  [XE]/dt,  and  performing  the  operation  [u]  d  [u]/dt  + 

[v]  d  [v]/dt.  Each  of  these  two  latter  terms  will  yield  two  eddy 
components  corresponding  to  the  last  two  terms  in  the  right-hand  side 
of  d  [X]/dt,  where  X  denotes  successively  [u]  and  [v].  Each  quadratic 
term  in  the  right-hand  side  of  Eqs.  5-i  and  5-j  yields  an  additional 
eddy  component,  so  that  the  six  eddy  components  in  (dXz/dt)*  are  as 
follows: 

=  pp  tan  0  [u*v*]  -  tan  0  [u*2]  - 

where  the  term  T  .  denotes  the  four  terms  in  the  individual  time 
e,3 

derivatives,  i.e., 


[u] 


jl 


9 


[u*v* ] 


m 


9 

"Sp 


[u*uf]>  +  [v] 


9 


[v*21 


m 


9 

T5 


[v*  urf  ]) 
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Noting  that  the  two  terms  involving  the  angular  momentum  [u*v* ]  can 
be  reduced  to  a  single  one,  i.e., 


¥{ 


tan  0  [u*v*]  -  m 


9  [u*v*  4=  -  iu,  2  _9_  [u*v*  1 


90  m 


and  using  9  [KE]/9t  =  0,  Eq.  B-6  reduces  to  the  form  of  Eq.  8-c  (Sec¬ 
tion  IX),  (Kurihara,  1970). 
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